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Tj- ; Abstract 

Based on the low-energy effective Hamiltonian with the generalized factor- 
\jQ ', ization, we calculate the new physics contributions to the branching ratios 

^ \ of the two-body charmless hadronic decays of B u and Bj mesons induced 

' by the new gluonic and electroweak charged-Higgs penguin diagrams in the 

general two-Higgs doublet models ( models I, II and III). Within the con- 
. sidered parameter space, we find that: (a) the new physics effects from new 

gluonic penguin diagrams strongly dominate over those from the new 7- and 
Q_i| Z°- penguin diagrams; (b) in models I and II, new physics contributions to 

O j' most studied B meson decay channels are rather small in size: from —15% to 

20%; (c) in model III, however, the new physics enhancements to the penguin- 
dominated decay modes can be significant, ~ (30 — 200)%, and therefore are 
• r- 1 , measurable in forthcoming high precision B experiments; (d) the new physics 

enhancements to ratios B{B — * Kr/') are significant in model III, ~ (35 — 70)%, 
and hence provide a simple and plausible new physics interpretation for the 
observed unexpectedly large B — ► Krf decay rates; (e) the theoretical predic- 
tions for B(B — ► K + tt) and B(B — > K°tt + ) in model III are still consistent 
with the data within 2a errors; (f) the significant new physics enhancements 
to the branching ratios of B -> if°7r°, K*rj, K* + ir-, K + <p, K*°lo, K* + <f) and 
K*°cj) decays are helpful to improve the agreement between the data and the 
theoretical predictions; (g) the theoretical predictions oiB(B — > PP, PV, VV) 
in the 2HDM's are generally consistent with experimental measurements and 
upper limits (90%C.L.) 

PACS numbers: 13.25.Hw, 12.15.Ji, 12.38.Bx, 12.60.Fr 
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I. INTRODUCTION 



The main objective of B experiments is to explore in detail the physics of CP violation, to 
determine many of the flavor parameters of the standard model (SM) at high precision, and 
to probe for possible effects of new physics beyond the SM Jl]-^]. Precision measurements 
of B meson system can provide an insight into very high energy scales via the indirect loop 
effects of new physics. The B system therefore offers a complementary probe to the searches 
for new physics at Tevatron, LHC and NLC colliders 

In B experiments, new physics beyond the SM may manifest itself, for example, in the 
following two ways |]]||: (a) decays which are expected to be rare in the SM are found to 
have large branching ratios; (b) CP-violating asymmetries which are expected to vanish or 
be very small in the SM are found to be significantly large or with a very different pattern 
with what predicted in the SM. These potential deviations may be induced by the virtual 
effects of new physics through loop diagrams. 

It is well known that the two-body charmless hadronic decays B — > h±h2 ( where hi and 
h 2 are the light pseudo-scalar (P) and/or vector(V) mesons ) play a very important role 
in studying CP violation and the heavy flavor physics 0,0. Several groups [§-|9j recently 
presented their systematic calculations for these B decay channels in the SM by using the low 
energy effective Hamiltonian [|T~0| P~^|| with the generalized factorization approach |7],|T^-T3 



Theoretically, the effective Hamiltonian is our basic tool to calculate the branching ratios 
and CP-violating asymmetry Aqp of B meson decays. The short and long distance quantum 
chromo-dynamical (QCD) effects in the hadronic decays are separated by means of the 
operator product expansion [lf|. The short- distance QCD corrected Lagrangian at next-to- 
leading order (NLO) is available now, but we still do not know how to calculate hadronic 
matrix element from the first principles. One conventionally resort to the factorization 
ansatz However, we also know that the non-factorizable contribution really exists 

and can not be neglected numerically for most hadronic B decay channels. To remedy 
factorization hypothesis, some authors |7| JT^fl5[1 introduced a phenomenological parameter 
N e ff (i.e. the effective number of color) to model the non-factorizable contribution to 
hadronic matrix element, which is commonly called the generalized factorization. On the 



other hand, as pointed by Buras and Silvestrini [17], such generalization suffered from the 



problems of gauge and infrared dependence since the constant matrix fy appeared in the 
expressions of effective Wilson coefficients C\ depends on both the gauge chosen and the 



external momenta. Very recently, Cheng et al. JT8[ studied and resolved above controversies 
on the gauge dependence and infrared singularity of C t - by using the perturbative QCD 
factorization theorem. In addition to the generalized factorization approach, a new approach, 
called as the QCD factorization []B5|, appeared recently JTP|pi]| , in which the decay amplitude 
is described by a kernel containing the 'hard' interaction given by a perturbatively evaluated 
effective Hamiltonian folded with form factors, decay constants and light-cone distributions 
of mesons into which the long distance effects are lumped. And some two-body hadronic 
B meson decays, such as B — > ttit and Kit modes, have been calculated in this approach 

EM- 

On the experimental side, CLEO collaboration reported the observations of thirteen 



B — > PP, PV decay channels and set new upper limits for many other decay modes [22-25 



The BaBar and Belle collaboration at SLAC and KEK also presented their first observation 
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for some P — > PP, PF decays at the ICHEP 2000 conference f2q , |27j| . Except for the decay 
channels P — > i^ry', the measured branching ratios for P — > /ii/i2 decays are generally in 
good agreement with the SM theoretical predictions based on the effective Hamiltonian 
with factorization. Unexpectedly large P — > Krj rate was firstly reported by CLEO in 1997 
p8| , and confirmed very recently by CLEO and BaBar Collaborations p3| , |29| , p6| . Although 
many possible mechanisms such as gluon and/or charm content in n' and the hairpin diagram 
have been considered in order to increase the theoretical predictions of B(B — > n'), it is now 
still difficult to explain the observed large rate for P — > Kv( decays [fj^,^,^]. This fact 
strongly suggests the requirement for additional contributions unique to the rj meson in the 
framework of the SM, or large enhancements from new physics beyond the SM. 

According to the studies in Refs. p0|-|3^| , we know that (a) an enhanced b — > sg can lead 
to a large B[B — > r]'X s ), and (b) the possible contributions to the ratio b — > sg in both type 
I and II two-Higgs-doublet models (2HDM) are not large enough to meet the requirement 
PH|3T| . Very recently, we calculated P^| , |33| the new physics enhancements to the branching 
ratios Bib — > sg) and B(b — > q'qq) with q' G {d, s} and q G {u, d, s} induced by charged- 
Higgs gluonic penguin diagrams in model III (the third type of 2HDM ) with inclusion of 
NLO QCD corrections p6| |, and found that the rate of b — > sg in model III can be enhanced 
significantly. The predicted charm multiplicity n c consequently become consistent with the 
measured n c , while the agreement between the theoretical predictions and the data of B$l 
is also improved by inclusion of the new physics effects. 

In this paper we calculate the new physics contributions to the branching ratios of ex- 
clusive two-body charmless hadronic decays P — > PP, PV, VV [] from new gluonic and elec- 
troweak charged-Higgs penguin diagrams in the general two-Higgs-doublet models (models 
I, II and III). We try to check the size and pattern of new physics effects on the exclusive 
two-body charmless B meson decays and to see if the new physics contributions in the model 
III can be large enough to provide the required enhancements for P — > Krf decay modes. 
We will present our systematic calculation of branching ratios for seventy six B — > /11/12 de- 
cay modes by employing the effective Hamiltonian with the generalized factorization 
We evaluate analytically all new strong and electroweak penguin diagrams induced by ex- 
changes of charged Higgs bosons in the quark level processes b — > qV* with q G {d, s} and 
V G {gluon, 7, Z}, and then combine the new physics contributions with their SM counter- 
parts and finally calculate the branching ratios for all seventy six exclusive P — > h\h2 decay 
modes. 

This paper is organized as follows. In Sec.[TT], we describe the basic structures of the 
2HDM's and examine the allowed parameter space of the general 2HDM's from currently 
available data. In Sec. |TTT| , we evaluate analytically the new penguin diagrams, combine 
the new physics contributions with their SM counterparts and find the effective Wilson 
coefficients . In Sec. W\, we present the formulae needed to calculate the branching 
ratios B(B — > h\h<i). In the following three sections, we calculate and show numerical results 
of branching ratios for P — ► PP, PV, and VV decay modes, respectively. We concentrate on 
those decay modes with well-measured branching ratios and sizable yields. The conclusions 



In the following, B always means B u or mesons. We here do not consider the decays of B s 
meson. 
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and discussions are included in the final section. 



II. THE GENERAL 2HDM AND EXPERIMENTAL CONSTRAINTS 



The simplest extension of the SM is the so-called two-Higgs-doublet models ||37|| . In 
such models, the tree level flavor changing neutral currents (FCNC's) are absent if one 
introduces an ad hoc discrete symmetry to constrain the 2HDM scalar potential and Yukawa 
Lagrangian. Lets consider a Yukawa Lagrangian of the form 

£y = vfjQi,L<PiU j>R + r)SQi,L<f>iD j!R + ^jQ.joJ'j.r + ^Q.j.o-iDjM + h.c, (1) 

where 0« (i = 1,2) are the two Higgs doublets of a two-Higgs-doublet model, 0i,2 = 2T20i 2 , 
Qi,L (Uj y R) with i = (1,2,3) are the left-handed isodoublet quarks (right-handed up-type 
quarks), Dj R are the right-handed isosinglet down-type quarks, while j^J and (i,j = 
1,2,3 are family index ) are generally the non-diagonal matrices of the Yukawa coupling. 
By imposing the discrete symmetry 



faD^-DM^TUi (2) 



one obtains the so called models I and II. In model I the third and fourth term in Eq.(|T]) 
will be dropped by the discrete symmetry, therefore, both the up- and down-type quarks 
get mass from Yukawa couplings to the same Higgs doublet 0i, while the 02 has no Yukawa 
couplings to the quarks. For model II, on the other hand, the first and fourth term in Eq.([l]) 
will be dropped by imposing the discrete symmetry. Model II has, consequently the up- and 
down-type quarks getting mass from Yukawa couplings to two different scalar doublets 0i 
and 2 . 

During past years, models I and II have been studied extensively in literature and tested 
experimentally, and model II has been very popular since it is the building block of the 
minimal supersymmetric standard model. In this paper, we focus on the third type of the 
two-Higgs-doublet model usually known as model III p8| , |39| . In model III, no discrete 



symmetry is imposed and both up- and down- type quarks then may have diagonal and/or 



flavor changing couplings with 0i and 02. As described in Ref. ||38|| , one can choose a suitable 
basis (H° , H 1 , H 2 , H^) to express two Higgs doublets |38j 



±( v/2x+ \ .1 / V2H+ 



(3) 



and take their vacuum expectation values as the form 

i> = 0, (4) 




where v = (^/2Gp) x l 2 = 246 GeV. The transformation relation between (H°, H 1 , H 2 ) 
and the mass eig enstates (IT, h°, A ) can be found in Ref. |38|. The H are the physical 



charged Higgs boson, H° and h° are the physical CP-even neutral Higgs boson and the A 
is the physical CP-odd neutral Higgs boson. After the rotation of quark fields, the Yukawa 



Lagrangian of quarks are of the form [38 
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£y 7 = Vi j Qi,L$iU jtR + n!'Quo.i)jM + &Qi,L$2Uj, R + BQ^faD^ + h.c, 



(5) 



where r^- D correspond to the diagonal mass matrices of up- and down-type quarks, while 
the neutral and charged flavor changing couplings will be |38|] 



fU,D 



V' 



-A 



^neutral 



u 

charged 



CKM, 



t 



V 



i ^charged 



VckmZ 



D 



(6) 



'1,2,3) are 



are free parameters to be determined by 



where Vckm is the Cabibbo-Kobayashi-Maskawa mixing matrix [40 
the generation index. The coupling constants Ajj 
experiments, and they may also be complex. 

In model II and setting tan/3 = 1)2/1)1 > 1 ( v\ and t>2 are the vacuum expectation values 
of the Higgs doublet fa and fa), the constraint on the mass M#+ due to CLEO data of 
b -> S7 [|T|] is M H + > 200 GeV for the charged Higgs boson in the 2HDM at the NLO level 
42fl . For model I, however, the limit can be much weaker due to the possible destructive 
interference with the SM amplitude. 

For model III, the situation is not as clear as model II because there are more free 
parameters here. As pointed in Ref. ||38|1 , the data of K° — K° and B® — B^ mixing processes 
put severe constraints on the FC couplings involving the first generation of quarks. Imposing 
the limit Ay = for j = (1, 2, 3) and assuming all other Ay parameters are of order 1, Atwood 



et al. (43| found a very strong constraint of M H + > 600 GeV by using the CLEO data of 
b — > 57 decay available in 1995. But this constraint can be lowered to M#+ > 400 GeV by 
using the new CLEO data of b — > 57 decay |35 



In Ref. [H3I , Aliev et al. studied the b 



S7 

decay in model III by extending the NLO results of model II [[42] to the case of model III, 
and found the constraint on the FC couplings. 

In a recent paper Chao et al. studied the decay b — > 37 by assuming that only 
the couplings X u = \\u\e ldt and A&& = | Afo^le*^ are nonzero. They found that the constraint 
on Mh+ imposed by the CLEO data of b — > sj can be greatly relaxed by considering the 
phase effects of X a and A^. From the studies of Refs. p5j|45| , we know that for model III 
the parameter space 



I A tt I 



0, for ij ^ tt, 
= 0.3, \Xbb\ = 35, 



or bb, 

e = (o° 



30 c 



M H+ = (200 ± 100) GeV, 



(7) 



are allowed by the available data, where 9 = 9^ — 9 t . 

From the LEP and Tevatron searches for charged Higgs bosons |H),17|, the new combined 
constraint in the (M#+,tan/3) plane has been given, for example, in Ref. [ I8| : the direct 



lower limit is M H + > 77GeV, while 0.5 < tan/3 < 60 for a relatively light charged Higgs 
boson with Mr+ ~ 100 GeV. Combining the direct and indirect limits together, we here 
conservatively consider the range of lOOGeV < M#+ < 300 GeV, while take = 200 

GeV as the typical value for models I, II and III. For models I and II we consider the range 
of 1 < tan (3 < 50, while take tan (3 = 2 as the typical value. In the following sections, 



We make the same ansatz on the & D couplings as the Ref. |38|. For more details about the 



definition of £ U,D one can see Ref. [38| 
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we calculate the new physics contributions to the exclusive two-body charmless decays of B 



meson in the Chao-Cheung-Keung (CCK) scenario of model III [45]. Model III in the CCK 
scenario has the following advantages: 

1. Since we keep only the couplings X a and X bb nonzero, the neutral Higgs bosons do not 
contribute at tree level or one-loop level. The new contributions therefore come only 
from the charged Higgs penguin diagrams with the heavy internal top quark. 

2. The new operators Og io and all flipped chirality partners of operators 0\ ... io as defined 



in Ref. |44[] do not contribute to the decay b — > 57 and the exclusive two-body charmless 
hadronic B decays under study in this paper. 

3. The free parameters are greatly reduced to \ u , X bb and M H + in model III, and tan/3 
and Mh+ in models I and II. 

III. EFFECTIVE WILSON COEFFICIENTS IN THE SM AND 2HDM'S 

In this section we evaluate the new gluonic and electroweak penguin diagrams and present 
the well-known effective Hamiltonian for the two-body charmless decays B — > hxhz with the 
inclusion of new physics contributions. For more details about the effective Hamiltonian 
with generalized factorization for B decays one can see, for example, Refs. J7H§. 

A. Operators and Wilson coefficients 

The standard theoretical frame to calculate the inclusive three-body decays b — > sqq 3 is 
based on the effective Hamiltonian [12,71 . 



G f I 



n eff (AB = 1) = -| j^c, (v ub v: s Q] + v cb v; s Q<) - v tb v* 



10 



CjQj + CgQg 

3=3 



where Cj and C g are Wilson coefficients, and the operator basis reads: 

Qi = {sq)v-A{qb)v-A, Q2 = {s a qp)v-A{qpb a ) V - A , 
with q = u and q = c, and 

Q3 = (sb) v ^ A ^2(q'q')v-A, Qa = (s a b ) v ^ A J2(l' Mv-a, 

q> q> 

Q5 = (sb) v ^ A J2(q'q')v+A, Qs = (s a b ) v ^ A J2(l' 'p&v+A, 

q< q> 

3 3 
Qr = -^{sb) v _ A ^e ql {q'q)v+ A , Qs = -{s a bp) v „ A J2 e q >(q' pq'Jv+A, 

q 1 q> 



(9) 
(10) 

(11) 

(12) 



3 For b — ► dqq decays, one simply makes the replacement s — > d. 
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3-3 
Q9 = ^(sb) v _ A J2e q/ (q'q') V -A, Qw = -z (s a bp)v-A e q'(<l' p<la)v-A, (13) 

q' q' 

Qg = ^m b s a a^(l + ^T^bpG^ (14) 

where a and (3 are the SU(3) color indices, T£g ( a — 1, ...,8) are the Gell-Mann matrices. 
The sum over q' runs over the quark fields that are active at the scale // = 0(m&), i.e., 
q' G {w, d, s, c, 6}. Qi and Q2 are current-current operators, C/3,4,5,6 and (57,8,9,10 are QCD 
and electroweak penguin operators, and Q g is the chromo-magnetic dipole (CMD) operator. 
Following Ref. 0, we also neglect the effects of the electro-magnetic penguin operator Q 77 , 
and do not consider the effects of the weak annihilation and exchange diagrams. 

In the SM, the Wilson coefficients C X (M W ), • • • , C 10 (M W ) at NLO level and C g (M w ) at 
leading order (LO) have been defined, for example, in Refs. [p| , p^| . The explicit expressions 
of the coefficients in the naive dimensional regularization (NDR) scheme can also be found 



easily in Refs. [11,12 



B. Contributions of the charged-Higgs penguin diagrams 



For the charmless hadronic decays of B meson under consideration, the new physics will 
manifest itself by modifying the corresponding Inami-Lim functions |^9| Cq(x), Dq(x), E q (x) 
and E' (x) which determine the coefficients C 3 (M W ), . . . , C W (M W ) and C g (M w ) in the SM. 
These modifications, in turn, will change for example the standard model predictions for the 
branching ratios of decays B — > hih 2 . The new strong and electroweak penguin diagrams can 
be obtained from the corresponding penguin diagrams in the SM by replacing the internal 
W lines with the charged-Higgs H lines, as shown in Fig.|l[ In the analytical calculations of 
those penguin diagrams, we use the dimensional regularization to regulate all the ultraviolet 
divergence in the virtual loop corrections and adopt the MS renormalization scheme. It is 
easy to show that all the ultraviolet divergence is canceled after summing up all Feynman 
diagrams. 

By evaluating analytically the new Z -, 7- and gluonic penguin diagrams induced by the 
exchanges of charged-Higgs boson H ± in the model III, we find the new Co, D , E and E' 
functions 



°0 



D III 



W III 



-x t 



16 



Vt 



l-Vt 



+ 



lit 



Vt)' 



■ IA 



tt\ 5 



-\H{y t )\\u\ 2 , 

-\i{yt)M , 



6 



J{y t )\\ tt \ 2 -Kiy^XuX^' 



(15) 

(16) 
(17) 
(18) 



with 



, 38y-79y 2 + 47y 3 Ay - 6y 2 + 3y\ 

H (y) = — ^ttti — — — + —nrr, — Hvl 



72(1 - yf 



12(1 - y y 



(19) 
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i(y) 
J(y) 

K(y) 



16y - 29y 2 + 7y 3 ^ 2y - 3y 2 

36(1 -yf 
2y + 5y 2 - y 3 

4(1-^ 
-3y + y 2 



6(1 -y) A 
3y 2 



2(1 -Z/) 4 



4(1 -y) 2 2(1-2/)' 



log[y] 



(20) 
(21) 
(22) 



where x t = m 2 /M^, ^ = m 2 /M] I+ , and the small terms proportional to m 2 /m 2 have been 
neglected. 

In models I and II, one can find the corresponding functions Co, D ,E and E' Q by 
evaluating the new strong and electroweak penguin diagrams in the same way as in model 
III 



n 1 — r 11 — Xt 
G °" G ° ~8tan 2 /3 



Ut 



+ 



Ut 



D> = D 



ii 



P / _ El// _ 

tan 2 /5 



3 tan 2 /3 



l-y t ■ (l-y t ) 2 
H(y t ) , 



1 



[J(j/ t ) - 6K(y t )\ 



E' 



ii 



3 tan 2 @ 



(23) 
(24) 
(25) 
(26) 
(27) 



where y t = m 2 /Mj I+ . 

We combine the SM part and the new physics part of the corresponding functions to 
define the functions at the scale /i = M\y as follows 



F (M ] 



w , 



+ ^0 



NP 



(28) 



where F e {Co, D , E , E' Q }. The explicit expressions of the functions Co, -Do, E and E' Q in 
the SM can be found, for example, in Ref. |i~2 |. 

Since the heavy new particles appeared in the 2HDM's have been integrated out at the 
scale M w , the QCD running of the Wilson coefficients Ci(M w ) down to the scale /x = 
O(mb) after including the new physics contributions will be the same as in the SM. By 



using QCD renormalization group equations [jTT] , |12|] , it is straightforward to run Wilson 
coefficients Ci(Mw) from the scale fi = 0(Mw) down to the lower scale \i = 0{rrib). Working 
consistently to the NLO precision, the Wilson coefficients C for i — 1, . . . , 10 are needed in 
NLO precision, while it is sufficient to use the leading logarithmic value for C g : 



C((x) = U(fi,M w )C(M ] 



C g (fi)= v u / 23 C g (M w )+Y / h iV a * 



(29) 
(30) 



where C(M W ) = (Ci(M w ), . . .,C 10 {M W )) T , U(fi,M w ) is the five-flavor 10 x 10 evolution 
matrix at NLO level as defined in Ref. pl]| , rj = a s (M^)/a s (fi), and the constants hi and 



a,- can also be found in Ref. [11 
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At the NLO level, the Wilson coefficients are usually renormalization scheme(RS) depen- 
dent. In the NDR scheme, by using the input parameters as given in Appendix and Eq.(^), 
and setting M H + = 200 GeV, 9 = 0° , tan/5 = 2 and fi = 2.5 GeV, we find the Wilson 
coefficients C e J* \\x) = C g + C 5 and Cj(/i) with i = 1, . . . , 10 in the SM and models I, II and 
III, and list them in Table |. From the numerical results as listed in Table |, one can easily 
see that: 

• The values of Ci(fi) ( i — 1, . . . , 10 ) in models I, II and III are almost identical with 
those in the SM. Only the coefficient C*'* in models II and III are clearly different 
from that in the SM. 

• It is the coefficient Cg'* partially induced by the new gluonic penguin diagrams which 
dominate the total new physics corrections to the decay processes under study. 



C. The effective Wilson coefficients 



We know that the unphysical RS dependence of Wilson coefficients will be cancelled by 
the corresponding dependence in the matrix elements of the operators in H-eff, as shown 



explicitly in Refs. |T2| , pO]| . Very recently, Cheng et al. |18| studied and resolved the so-called 
gauge and infrared problems fl7| of generalized factorization approach 0. They found that 
the gauge invariance is maintained under radiative corrections by working in the physical on- 
mass-shell scheme, while the infrared divergence in radiative corrections should be isolated 
using the dimensional regularization and the resultant infrared poles are absorbed into the 
universal meson wave functions 



IS 



The one-loop matrix elements can be rewritten in terms of the tree-level matrix elements 
of the effective operators [|7| 



( Sq iqi\H en \b) = Y.Ct ti ^)(s q >qi\0 3 \b) 



\ tree 



(31) 



'■j 



where C- ff (/j) (i = 1, . . . , 10) are the effective Wilson coefficients. In the NDR scheme and 
for SU(3)c, the effective Wilson coefficients C e ^ can be written as [0,0, 



C 



eff 



1 + ^ (r y + 7y log- 



J IJ 



Q:. s 



24vr 



a, 



4 (C t + C p + C g ) + B[C 



71 



(32) 



where A[ = (0, 0, -1, 3, -1, 3, 0, 0, 0, 0) T , B[ = (0, 0, 0, 0, 0, 0, 1, 0, 1, 0) T , the matrices r v and 
7v contain the process- independent contributions from the vertex diagrams. Like Ref. 0, 
we here include vertex corrections to C-j — Cio []. The anomalous dimension matrix 7y has 
been given explicitly, for example, in Eq.(2.17) of Ref. 0. Note that the correct value of 
the element (r^DRlm and (rNDFijss should be 17 instead of 1 as pointed out in Ref. [^TJ, fy 
in the NDR scheme takes the form 



4 The reliability of the generalized factorization approach is improved by this progress. 
5 Numerically, such corrections are negligibly small. 
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3 ) 



(33) 



1.2 



The function Ct 

Qs,.. 



C p , and C g describe the penguin-type corrections to the operators 
and the tree-level diagram of the operator Q g respectively. We here follow 

The effective Wilson coefficients Cf 



Q 

the procedure of Ref. fl5| to include C g in (^ 



m 



and infrared safe. The explicit expressions of functions C t , C p , and C g in the NDR scheme 



Eq.([32]) are now scheme and scale independent in NLO precision, and also gauge invariant 
and infrared safe. The explicit expressions 
have been given, for example, in Refs. |7].[J 

C t 



~ + ^G(m u ) + ^G(m c ) 



Ci, 



Cn 



a, 



- - G{m q ) - G(m b ) 



c 3 + E 



i=u,d,s,c,b 



3 - gk; 



(C 4 + c 6 



2m b r eS 



A, 



+ -^G(m c 



(Ci + 3C 2 ) 



with A r , 



V q >bV*, q . The function G(m) is of the form |5 



where z 



G{m) 
k 2 /(4m 2 ), and 
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--lnl 
3 1 



+ 



2/i 2 _ 2(1 + 2z) 
3m 2 3z 




- in 



z < 1, 

z > 1. 



(34) 
(35) 

(36) 
(37) 

(38) 
(39) 



where k is the momentum transferred by the virtual gluon, photon or Z to the q'q' quark pair 
in the inclusive three-body decays b — > qq'q', and m is the mass of internal up-type quark 
in the penguin diagrams. For k 2 > 4m 2 , an imaginary part of g(z) will appear because of 



the generation of a strong phase at the uu and cc threshold []52|-[54[ . 

For the two-body exclusive B meson decays any information on k 2 is lost in the fac- 
torization assumption, and it is not clear what "relevant" k 2 should be taken in numerical 
calculation. One usually uses the "physical" range for k 2 : m 2 /4 ~ k 2 ~ ml/2. Follow- 
ing Refs. we also use k 2 = ml/2 in the numerical calculation and will consider the 
fc 2 -dependence of branching ratios of charmless B meson decays for several typical decay 
channels. 



10 



IV. DECAY AMPLITUDES IN THE BSW MODEL 



In numerical calculations, two sets of form factors at the zero momentum transfer from 
the Bauer, Stech and Wirbel (BSW) model fll3| , as well as Lattice QCD and Light-cone QCD 



sum rules (LQQSR) |5l| will be used respectively. Explicit values of these form factors can 
be found in Ref. and have also been given in Appendix. Following Ref. JFj, the seventy 
six decay channels of B u and Bj mesons are classified into five classes according to their 
JV e / / — dep endence : 

• Class-I: including four decay modes, B° — > tt~tt + , p ± 7r T and B° — > p~K + , the large 
and N e H stable coefficient a\ plays the major role. 

• Class-II: including ten decay modes, for example B — > 7T 7T , and the relevant coeffi- 
cient for these decays is a 2 which shows a strong N e ^ -^-dependence. 

• Class-Ill: including nine decay modes involving the interference of class-I and class-II 
decays, such as the decays B + — > 7r + ?/. 

• Class-IV: including twenty two B — > h±h 2 decay modes such as B — > Krf> decays. The 
amplitudes of these decays involve one (or more) of the dominant penguin coefficients 
a 4,6,9 with constructive interference among them. The Class-IV decays are N e ^ stable. 

• Class- V: including twelve B — ► hih 2 decay modes, such as B — > 7r ?/') and B — > (pK de- 
cays. Since the amplitudes of these decays involve large and delicate cancellations due 
to interference between strong N e ^- dependent coefficients a^^j^Q and the dominant 
penguin coefficients a 4j6j 9, these decays are generally not stable against N e ^ . 



With the factorization ansatz [|T^,[56|,[)7]] , the three-hadron matrix elements or the decay 
amplitudes < XY\H e ff\B > can be factorized into a sum of products of two current matrix 
elements < X|Jf|0 > and < Y\ J 2fl \B > ( or < Y\j£\0 > and < X\ J 2fl \B >). The explicit 
expressions of the matrix elements in terms of decay constants {fxiQx) an d the Lorentz- 
scalar form factors A 0> i t2 (k 2 ) and F ^(k 2 ) can be found, for example, in Refs. pT5|J58| J?[] . 



In the B rest frame, the branching ratios of two-body B meson decays can be written as 
B(B -> XY) = r B ^r\M(B - XY)\ 2 (40) 
for B — > PP decays, and 



B(B -> XY) = r B ^\M(B - XY)/(e ■ p B )\ 2 (41) 



v 



for B — > PV decays. Here t(B~) = 1.65 ps and t(B%) = 1.56 ps [ 59| , p B is the four- 
momentum of the B meson, My and e is the mass and polarization vector of the produced 
light vector meson respectively, and \p\ is the magnitude of momentum of particle X and Y 
in the B rest frame 



H = wrJ {Ml ~ {Mx + My ^ m * - ( Mx - My ^ ■ ^ 
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For B — > VV decays, one needs to evaluate the helicity matrix elements H\ =< 
V\(A)V" 2 (A)|iJ e// | J B) > with A = 0,±1. The branching ratio of the decay B -> V X V 2 is 
given in terms of H\ by 

B(S -> V 1 V 2 ) = r B -^L [|^ |2 + \ H+X \2 + i^^l _ (43) 



8ttM 



The three independent helicity amplitudes H , H +1 and H_i can be expressed by three 
invariant amplitudes a, 6, c defined by the decomposition 



H x = ze"(A)^(A) 



ZC rv fl 



(44) 



where p\ }2 an d M 12 are the four momentum and masses of Vip, respectively, p = p\ + p 2 is 



the four-momentum of B meson, and 



H ±1 — a ± cVx 2 — 1, -ffo — — b (x 2 — 1) (45) 



Ml-M 2 -M 2 fl . 

x = — 46 

2M 1 M 2 K 1 

In the generalized factorization ansatz, the effective Wilson coefficients Cf will appear 
in the decay amplitudes in the combinations: 

Q e ff (J e ff 

fl2i-i = C 2i _ 1 + jt^jj, «2i = C 2i + —^jj , (i = 1, . . . , 5) (47) 

where the effective number of colors N e ^ is treated as a free parameter varying in the range 
of 2 < N e ff < oo, in order to model the non-factorizable contribution to the hadronic 
matrix elements. It is evident that the reliability of generalized factorization approach has 
been improved since the effective Wilson coefficients C\ * appeared in Eq. (f|7|) are now gauge 
invariant and infrared safe. Although N e ^ can in principle vary from channel to channel, but 
in the energetic two-body hadronic B meson decays, it is expected to be process insensitive as 
supported by the data [|. As argued in Ref. @, N eff (LL) induced by the (V — A)(V — A) 
operators can be rather different from N e ^(LR) generated by (V — A)(V + A) operators. In 
this paper, however, we will simply assume that N e ^(LL) = N e f* (LR) = N e ^ and consider 
the variation of N e ^ in the range of 2 < N e ^ < oo since we here focus on the calculation of 
new physics effects on the studied B meson decays induced by the new penguin diagrams in 
the two-Higgs-doublet models. For more details about the cases of N e ^(LL) ^ N e ^(LR), 
one can see for example Ref. 0. We here will also not consider the possible effects of final 
state interaction (FSI) and the contributions from annihilation channels although they may 
play a significant rule for some B meson decays. 

Using the input parameters as given in Appendix, and assuming k 2 = ml/2, M H + = 
200 GeV, 9 = 0° and tan/? = 2, the theoretical predictions of effective coefficients a« are 
calculated and displayed in Table and Table III for the transitions b — > d (b —>■ d ) and 



b — ► s (b — > s), respectively. For coefficients a 3 , . . . , aio, the first, second and third entries in 
tables (mplD refer to the values of in the SM and models II and III, respectively, a, in 
model I are very similar with those in the SM and hence was not given explicitly. 
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All branching ratios in the following three sections are the averages of the branching ratios 
of B and anti--B decays. The ratio 53 describes the magnitude of new physics corrections 
on the SM predictions of the decay ratios and is defined as 

6B(B - XY) - B{B ~* XY)2HDM ~ B{B ~* XY)SM (48) 
5B(B - XY) - xy)SM (48) 



V. -B — > PP DECAYS 

Using formulae as given in last section, it is straightforward to find the decay am- 
plitudes of B — > PP decays. As an example, we present here the decay amplitude 

M(B~ -> 7T-7T ) =< TT-TT°\H eff \B- >, 

G 



M{B- -> tt-tt ) = -f {v ub v: d ( ai M:j° + a 2 Mi:°) 



with 



2 



'a 4 + aio + (a 6 + a 8 )i?i) Ml * 



«4 + ^(a 7 - « 9 ) - ^ + (a 6 - |)i* 2 ) M^/ 



(49) 



i?2 



Ml/ 



2mf 



(m b - m u ){m u + m rf ) 



m d (m b - m d ) 



-zlm 



B 



-Km 



B 



(mlo) 



(50) 

(51) 

(52) 
(53) 



where f n is the decay constant of 7r meson. The form factor F^^(m 2 ) can be found in 
Appendix. Under the approximations of setting m u = rrid and m n o = m n -, the decay 
amplitude M(B~ — ► tt^tt ) in Eq. (^9|) will be reduced to the form as given in Eq.(80) of 
Ref. 0]. In the following numerical calculations, we use the decay amplitudes as given 
in Appendix A of Ref. directly without further discussions about details of individual 
amplitude. 

In tables |TV| - |VI| , we present the numerical results of the branching ratios for the twenty 
B — > PP decays in the framework of the SM and models I, II and III by using the BSW and 
LQQSR form factors, respectively. Theoretical predictions are made by using the central 
values of input parameters as given in Eq.(|7|) and Appendix, and assuming M H + = 200 
GeV, 6 = 0°, tan/5 = 2 and N e ^ = 2, 3, oo in the generalized factorization approach. The 
/c 2 -dependence of the branching ratios is small in the range of k 2 = m 2 /2 ± 2 GeV 2 and 
hence the numerical results are given by fixing k 2 = ml/2. 

The currently available CLEO data f2^24| are listed in the last column of table uM. 
From the numerical results, we see that: 
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For B — > Kr]' decays, the observed branching ratios are clearly much larger than the 



SM predictions [g5|,g]|j. All other estimated branching ratios in Table [IV] are, however, 



consistent with the new CLEO, BaBar and Belle measurements or upper limits. 

In model III, the new physics corrections to most class-II, IV and V decay channels 
can be rather large and insensitive to the variations of the mass M#+ and the color 
number N eff : from 20% to 90% w.r.t the SM predictions for both cases of 6 = 0°, 30°. 
For tree-dominated decay modes B — > n + ix~ , 7r + 7r°, n + r]('\ the new physics corrections 
are small in size. 

In models I and II, however, the new physics corrections to all B — > PP decay modes 
are small in size within the considered parameter space: less than 3% in model I, and 



~ (—20 — 0)% in model II, as shown in Table [Vl|. So small corrections will be masked 
by other large theoretical uncertainties. 

In model III, the new gluonic penguins will contribute effectively through the mixing 
of chromo-magnetic operator Q g with QCD penguin operators Qs — Qq. The C e J^ will 
strongly dominate the new physics contributions to all B — * h\\i2 decay modes. 

The central values of the branching ratios obtained by using the LQQSR form factors 
will be increased by about 15% when compared with the results using the BSW form 
factors, as can be seen from Table |V|. We therefore use the BSW form factors only to 
calculate the new physics effects on the ratios BiB — > hih 2 ) and treat the difference 
induced by using different set of form factors as one kind of theoretical uncertainties. 



A. B — > 7T7T, Kit decays 

There are so far seven measured branching ratios of B — > PP decays: one B — > 7r + 7T 
decay, four B — > Kir and two B — > Kr]' decays ]2"B|,2"3,|2"E, 



BiB - Tr+vr-) - \ (43 ^ ± °- 5) X 10_6 [CLE0] ' (54) 

7T ) - I (9 3+2.8+1.2) x 1Q -6 [BaBar]; m 

BiB - K+rr ) - I ( 1L6 ^^) >< 10 " 6 t CLE °]' (ss) 

vr ) - j (lg 8 +5.5 ± 2 3) x 10 _ 6 [Belle]> (55) 

f (17.2±fj ± 1.2) x 10~ 6 [CLEO], 

B{B -> K + w-) = < (12.5±|:g ± 2.3) x 10~ 6 [BaBar], (56) 
[ (17A±H ±3.4) x 10- 6 [Belle], 

B{B -> K\ + ) = (18.2^o ± 1.6) x lO" 6 [CLEO], (57) 

BiB - K%°) - I (1L(i X 10_6 [ CLE °]' (58) 

/ — j (21±? : |±i:l) x 10- 6 [Belle], (58j 

Bi B - K + n') - I (8 °^° ± 7) X 10_6 [CLE ° ] ' (59) 

V } ~ \ (62 ± 18 ± 8) x 10- 6 [BaBar], (59j 

B{B -> K%') = (89±il ± 9) x 10~ 6 [CLEO]. (60) 
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The measurements of CLEO, BaBar and Belle Collaborations are in good agreement 
with each other within errors. These decays are sensitive to the relevant form factors 
F^, F^ n , F^ v \ etc., and to the value of N e ^ '. 

As a Class-I decay channel, the B° — > tt + tt~ decay is dominated by the b — > u tree 
diagram. The band between two dots lines in Fig|| shows the CLEO measurement. Since the 
new physics corrections are very small in size, less than 3% within the considered parameter 
space, the four curves for the SM and 2HDM's are close together and can not be separated 
clearly. The theoretical predictions look higher than the CLEO measurement, but they are 
still consistent with BaBar measurement because of very large error of BaBar data. In fact 
the theoretical predictions for B(B — > 7r + 7r _ ) in the SM and 2HDM's are still consistent 
with the CLEO data at the 2a level if we consider currently still large theoretical and 
experimental uncertainties. On the other hand, if we take the average of CLEO and BaBar 
measurements, B(B® — > 7r + 7r~) = (5.5 ± 1.5) x 10~ 6 , as the experimental result, then the 
constraint on F<f ^(0) from the data will be F^ n (0) = 0.25 ± 0.03 by setting A = 0.2205, 
A = 0.81; p = 0.12, r] = 0.34, N e ^ = 3, and by neglecting FSI also. 

In the SM, the four Class-IV decays B — ► Krr are dominated by the b —>■ sg gluonic 
penguin diagrams, with additional contributions from b — > u tree and electroweak penguin 
diagrams. Measurements of B — > Ktc decays are particularly important to measure the 
angle 7. In model III, the new physics enhancements to the branching ratios B(B — > Kit) 
are significant, ~ (50 — 60)%, and show a moderate dependence on the variations of other 
parameters, as illustrated in figures (0-0). In models I and II, however, the new physics 
corrections are always very small in size. 

For the decays B — > K + 7t~ , K°it + , the theoretical predictions in model III are higher 
than the CLEO data as shown in Fig.[| and Fig.|5|, but they are still consistent with the 
CLEO data at the 2a level if we consider currently still large theoretical uncertainties. As 
a simple illustration of effects of the theoretical uncertainties, we recalculate the branching 
ratios of B — > K + n~ and K°ir + decays by using F B ^ 7r (0) = 0.25 instead of the ordinary 
BSW value F^^fo) = 0.33 while keeping all other input parameters unchanged, and find 
numerically that 

7T j - I (1? 2+4.3+3.7) x 10 -6 in Modd ni) W 

B(B - - { X 10 ~ 6 in SM ' f62) 

I (19.9±|;°±li) x 10- 6 in Model III, {b2} 

for N e -ff = 3 and M H + = 200 GeV. Here the first and second error correspond to F^ n (0) = 
0.25 ± 0.03 and 2 < N e ^ < 00, respectively. It is evident that the theoretical predictions of 
the two ratios B(B -> K + tt~) and B(B -> K°tc + ) in the SM and model III can lie within 
the CLEO data. 

Figures (§-(§) show the mass and iV e ^-dependence of the branching ratios for four B — > 
Kit decay modes in the SM and models I, II and III, using the input parameters as given 
in Eq.(0) and Appendix, and assuming 6 = 0°, tan/5 = 2, and k 2 = ml/2. For Figs. da- 
la), we set iV e// = 3 and assume that M H+ = 100 - 300GeV. For Figs.(|b-|b), we set 
Mh+ = 200GeV, and assume that 1/N e ^ = — 0.5. In all four figures, the band between 
two dots lines shows the corresponding CLEO measurements with 2a errors. For B — > K°7r° 
decay, the inclusion of new physics contribution will improve the agreement between the 
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data and theoretical prediction, as illustrated in Fig.^. For other three B — > decays, the 
theoretical predictions in the model III are still consistent with the data if the theoretical 
uncertainties are taken into account. 



B. B — > Krf> decays and the new physics effects 

For B + — > K + f) and B° — > K°rj decay modes, the new physics corrections are large 
(small) in model III (models I and II). The theoretical predictions in the SM and 2HDM's 
are consistent with the new CLEO upper limits. 

For B — > Kvj decay modes, the situation is very interesting now. In 1997, CLEO firstly 



reported the unexpectedly large B — > Krf rates [28|, which is confirmed very recently by 
CLEO with the full CLEO II/II.V data sample of 19 million produced B mesons @|29[. 
The Krf signal is large, stable and has small error (~ 14%). Those measured ratios are 



clearly much larger than the SM predictions as given in table ([TV|). In [51], Cheng and Yang 
considered various possible enhancements to Kr( decay modes in the framework of the SM 
0, but found that the net enhancement is not very large: B(B ± — > K rj') = (40 — 50) x 10~ 6 , 
which is smaller than the CLEO data [] At present, it is indeed difficult to explain the 
observed large rate for B — > Krf This fact strongly suggests the requirement for 

additional contributions unique to the rf meson in the framework of the SM, or from new 
physics beyond the SM. 

In models I and II, the new physics contributions are too small (or negative) to provide 
the required enhancement. This feature remains unchanged within the considered range of 
tan/5 = 1 - 50. 

In model III, however, the new physics enhancements are significant, ~ 60%, and have 
a moderate dependence on M H + and N e ^ , as illustrated by the solid curves in Figs. (J7|-[rDD 
where only the central values of theoretical predictions in model III are shown. If we take into 
account other theoretical uncertainties, the theoretical predictions for ratios B{B — > Kr/ 1 ) 
in model III will become consistent with the CLEO data: 



(63) 

B(B° -> K°rj') = { (19 - 53) x 10~ 6 [SM] , (64) 



(69 


- 92) x 10~ 6 


[CLEO] , 


(20 


- 52) x 10~ 6 


[SM] , 


(34 


- 74) x 10~ 6 


[Model III] 


(71 


- 109) x 10- 


5 [CLEO] , 


(19 


-53) x 10~ 6 


[SM] , 


(33 


- 73) x 10- 6 


[Model III] 



6 As discussed in ||, B — > Krj decay may get enhanced due to (i) a small m s at the scale m&, 
(ii) the sizable SU(3) breaking, (iii) large Fq^ 11 , (iv) the rf charm content, and (v) constructive 
interference in tree amplitudes. But these possible enhancements are partially washed out by the 
anomaly effect in the matrix element of pseudoscalar densities [ p^j3^| . 

7 Although this prediction is consistent with BaBar measurement, one should note that the error 
of BaBar measurement is still much larger than that of CLEO data. More statistics is clearly 
required for BaBar to make a definite conclusion. 
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Here the major theoretical uncertainties induced by using different set of form factors and 
varying k 2 , rj and N e ^ in the ranges of 5k 2 = ±2GeV 2 , Srj = ±0.08 and N e ^ = 2 — oo have 
been taken into account. 

Figures ([|,[|) show the mass and N e ^ dependence of BiB — > Krj') in the SM and 
2HDM's. The upper dots band shows the CLEO measurements with 2a errors. The short- 
dashed, dot-dashed, long-dashed and solid curve refers to the theoretical predictions in 
the SM, models I, II and III, respectively. As shown explicitly in Figs.(|S|,|T(lD, in which 
the short-dashed, long-dashed and solid curve corresponds to the model III predictions for 
jye// _ 2, 3, oo respectively, the theoretical predictions become now consistent with the 
CLEO measurement due to the inclusion of new physics enhancement in model III. 



VI. B — > PV DECAYS 

In tables ( |V11[ - |LXD we present the branching ratios for the thirty seven B — > PV decay 
modes involving b — > d and b — > s transitions in the SM and models I, II and III by using 
the BSW form factors and by employing generalized factorization approach. Theoretical 
predictions are made by using the same input parameters as those for the B — > PP decays 
in last section. 

For studied thirty seven B — > PV decays, two general features are as follows: 

• The theoretical predictions for those seven measured decay rates are consistent with 
the CLEO data within 2a errors. All other estimated branching ratios in the SM and 
2HDM's as given in tables (|yiiyiX| ) are all consistent with the new CLEO upper limits. 



• For most decay modes, the differences induced by using whether BSW or LQQSR 
form factors are small, ~ 15%. We therefore use the BSW form factors only in the 
calculation of new physics effects. 

There are so far seven measured branching ratios of B — > PV decays. For the first 
three decay modes, B — > p ± p =F , p°vr + , um + , the new physics corrections are small in size, 
< 5%, and have a weak dependence on M#+ and N e ^ , as shown in tables (|Vlll|jrX| ). 
Consequently, the theoretical predictions in the SM and models I, II and III agree well with 
CLEO measurements. 

Because of the appearance of very large B{B — > Kr) 1 ), the decay modes B — > K*v^ ~> also 
draw more attentions now. Very recently, CLEO and Belle reported their first observation 



2,23,27] of B ->• K*r], K* + 7r~ and B K + <p decays: 



B(B + -> K* + r]) = (26.4±|;| ± 3.3) x 10~ 6 , (65) 

B{B° -> K*° v ) = (13.8±1;1 ± 1-6) x 10" 6 , (66) 

B(B° -> K* + n-) = (22±|±|) x 10~ 6 , (67) 

B(B + -> K + 4>) = (17.2l|;| ± 1.8) x 10" 6 , (68) 



while the theoretical predictions in the SM and model III are 
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B(B+ 


-> K* + r]) 


B(B + 




B{B + - 




B(B + 


^if+0) 



(2 


-4) x iir 6 


[SM] , 


(2 


- 5) x 1(T 6 


[Model III] 


(2 


- 5) x 10~ 6 


[SM] , 


(3 


-6) x 10~ 6 


[Model III] 


(7 


- 16) x 10- 


6 [SM] , 



(69) 
(70) 

(10 - 22) x 10~ 6 [Model III] , ^ 
(0.5 - 28) x 10- 6 [SM] , 

(1-39) x 10~ 6 [Model III] , 1 ' 

where the uncertainties induced by using the BSW or LQQSE form factors, and setting 
k 2 = ml/2 ± 2GeV 2 , rj = 0.34 ± 0.08 and N e ^ = 2 — oo, have been taken into account. 
Although the central values of the theoretical predictions in the SM are much smaller than 
the corresponding central values of the CLEO measurements, the theoretical predictions 
are still consistent with the data within 2a errors because current experimental error is 
still large. Further improvement of experimental measurements about the decay modes 
B — > K*T], K* + tt~ will tell us whether there is any discrepancy between the theory and 
experiments for these three decay modes. At present, any positive contributions to the 
above three branching ratios from new mechanisms in the SM or from new physics beyond 
the SM are clearly preferred by the CLEO data. 

In models I and II, the new physics contributions are small in size: from —15% to 20% for 
most B — > PV decay modes, and have weak dependence on M#+, tan/3 and N e ^ , as shown 



in table (pC|), and illustrated in Figs.(p4]-|T4]) where the long-dashed line shows the theoretical 
predictions in the model II []. This feature remains unchanged within the considered range of 
tan f3 — 1 — 50. When tan j3 becomes larger, the size of new physics corrections will become 
even smaller. 

In model III, however, the new physics contributions are significant, from 30% to 110%, 
and have also weak dependence on M#+, 9 and N e ^ . These new physics enhancements are 
very helpful to improve the agreement between the theoretical predictions and the data, as 
shown in Eqs.( |69] - |72|) and illustrated in Figs. (|TT]-^4[) 



Figs.( |TT| - |T4] ) show the mass and N e ^ dependence of the branching ratios for B — > 
K* + i], K*°7], K* + tc~ and B — > K + <p decays. The dot-dashed line is the SM prediction, 
while the long-dashed and solid curve corresponds to the predictions in models II and III, 
respectively. The theoretical uncertainties are not shown in these figures. The dots band in 
Figs.( |Tll - |l3"D (Fig.[L4]) corresponds to the CLEO data with la (2 — a) error. 

From Fig. (|TTD , it can be seen that the CLEO measurement of the ratio B{B + — > K* + if) is 
much larger than theoretical predictions in the SM and 2HDM's. More positive contributions 
to this decay mode are needed to improve the agreement between the data and theoretical 
prediction. For B — > K*°r] decay, the inclusion of new physics contribution in the model 
III leads to a better agreement between data and theory if we take into account still large 
theoretical uncertainties. For B — > K* + n~ and K + (p decays, the theoretical prediction 



"Because the lines for the SM and model I are too close to be separated clearly, we do not draw 
the line for model I in all four figures for B — > PV decays. 
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becomes now consistent with the CLEO and Belle measurements within la error due to the 
large new physics enhancement in the model III. 

For B — > K* + i]' and K* t]' decays, the new physics contributions in model III are large in 
size, from —77% to 200%, as shown in Table [VII]] . But the theoretical predictions for these 
two decay modes are A e ^-dependent and still far below the current CLEO upper limits. 



VII. B VV DECAYS 



Using the formulae as given in Sec. [IV], it is straightforward to calculate the branching 
ratios of nineteen B u ^ — > VV decays. As an example, we show here the calculation of 
the branching ratio for the Class-V decay B~ — > p~u (b — > d transition). We firstly find 
the explicit expressions of the helicity amplitude H\ =< p~ (\)u(\)\H e ff\B~ >, and then 
compare this amplitude with the standard form as defined in Eq . ([44]) to extract out the 
process dependent coefficients a, b and c 



1 

V2 



fi f u M u (M B + M p )A^ p {Ml) + f 2 f p M p (M B + M ul )Af- ¥U {M i 



B- 



b = fi- 
c = fi- 



M B + M p 



M B + M^ 



>(M- 



I)' 



V2f^M PvB ^ {M2j + h ^Af M 



M B + M P 



M B + M u 



with 



G p 
72 



V ub V* d a 2 - V tb V* d (2a 3 + a 4 + 2a 5 + ^a 7 + ^a 9 - ^a w 
[VubV*^! - V tb V* d (o 4 + aio)] , 



(73) 
(74) 
(75) 

(76) 
(77) 



where the coefficients ai v .. j io have been defined in Eq . (|47|) , the form factors and other input 
parameters can be found in Appendix. With these coefficients a, b and c, the branching ratio 
B(B~ — > p~u>) can finally be written as 



B{B~ 



P 



\P\ 



B - 8vrM| 



o\ 2 + \H +1 \ 2 



+ l^-i| 2 ) 



(78) 



where \p\ and Hi^-i have been given in Eqs.(|42|) and (|45|) . 

In tables (JX|- |X-1J] ) we present the branching ratios for the nineteen B — > VV decay 
modes involving b — > d and b —> s transitions in the SM and models I, II and III. Theoretical 
predictions are made by using the same input parameters as those for the B — > PP, PV 
decays in last two sections. 

For B — > VV decay modes, the differences induced by using whether BSW or LQQSR 
form factors are around ten percent in the SM and models I, II and III. We therefore show 
the numerical results obtained by using the BSW form factors only for the cases of models 
I, II and III. For all nineteen B — > VV decays under study, the theoretical predictions in 
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the SM and 2HDM's are still under or far away from the current CLEO upper limits, as can 
be seen from tables (|X- XH ). 

In models I and II, the new physics contributions to B — > VV decays are small in size: 
from —15% to ~ 10% as shown in Table [XII], and therefore will be masked by other large 
theoretical uncertainties. This feature remains unchanged within the considered range of 
tan /3 — 1 — 50. When tan/3 becomes larger, the size of new physics corrections will become 
smaller. 

In model III, however, the new physics contributions to different channels are varying 
greatly: from -11% to ~ 110%, assuming M H + = 200GeV, N e " = 2 - oo and 9 = 0° - 30°. 
For decay modes B — > K*°u, K* + <ft, K*°(f>, for example, the new physics enhancements are 
significant: ~ (60 — 110)%. And hence the theoretical predictions in model III are close 
to or slightly surpass the current CLEO upper limits, as illustrated in Figs. ( |15l - |Hj ) where 
the upper dots line shows the corresponding CLEO upper limits at 90%C.L.. These decay 
modes will be observed soon. 

Figures (|T5| - |T7|) show the mass and N e ^ dependence of the ratios BiB — > K*°uj) and 
B(B — > K*(f). The dot-dashed line is the SM prediction, while the long-dashed and solid 
curve corresponds to the predictions in models II and III, respectively. As the Class-V 
decays, these three decays show strong N e ^ dependence as illustrated in Figs. ( |T5| - P~7[ ) . 



VIII. SUMMARY AND DISCUSSIONS 

In this paper, we calculated the branching ratios of two-body charmless hadronic B 
meson decays B u ^ — > PP, PV, VV in the SM and the general two-Higgs-doublet models by 
employing the NLO effective Hamiltonian with the generalized factorization. 

In Sec. ([1], with the help of previous works p8| , ^3| - ^ , p5| 1 , we gave a brief review about 



the 2HDM's and studied corresponding experimental constraints on models I, II an III. In 



Sec. |Tl|, we evaluated analytically all new gluonic and electroweak charged-Higgs penguin 
diagrams and found the effective Wilson coefficients Cf in the SM and models I, II and 
III. In Sec. |V|, we presented the formulae needed to calculate the branching ratios BiB — > 
PP,PV,VV). 



In sections |V|- |VII| , we calculated the branching ratios for seventy six B — > PP, PV, VV 



decays in the SM and models I, II and III, presented the numerical results in tables (|TV| - |XHD 



and displayed the Mh+ and Af e -^-dependence for several phenomenologically interesting 
decay modes in Figs.(|2]-|T7|). 

From the numerical results, we find following general features about the new physics 
effects on the exclusive charmless hadronic B — > PP, PV, VV decays studied in this paper: 

1 The SM predictions for the B meson decay rates presented in this paper agree well 
with those appeared in Refs. |[7]||. 

2 The new physics effects due to new gluonic penguin diagrams strongly dominate over 
those from the 7— and Z°— penguin diagrams induced by exchanges of charged-Higgs 
bosons appeared in models I, II and III. 
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For models I and II, the new physics contributions to the decay rates BiB — > hxh-i) are 
always small in size: from —15% to 20% for most decay modes. So small contributions 
will be masked by other still large theoretical uncertainties. 

For model III, however, the new physics enhancements to penguin-dominated decay 
modes can be significant, ~ (30 — 200)%, and therefore can be measured in high 
precision B experiments. In general, the new physics contributions in model III are 
large (small) for penguin- dominated (tree-dominated) B meson decay channels. 

The uncertainties of the theoretical predictions for the branching ratios of B — > hih^ 
decays induced by varying k 2 , r], 9, tan/5, and M H + are varying from ~ 10% to ~ 50% 
within the range of k 2 = m 2 b /2 ± 2GeV 2 , rj = 0.34 ± 0.08, 9 = 0° - 30°, tan/3 = 1 - 50, 
and M H + = 200 ±100 GeV. The dependence of decay rates on whether using the BSW 
or LQSSR form factors are weak, ~ 10%. The N e ^ — dependence of branching ratios, 
however, are varying greatly for different decay modes. 

For phenomenologically interesting B — > Krf decay modes, the new physics enhance- 
ments are significant in model III: ~ (35 — 70)%, and have a moderate dependence 
on M H + and N e ^ . The theoretical predictions for B(B — > Kr}') therefore turn to be 
consistent with the CLEO data in model III, as illustrated in Figs.(|8|,|l(]). For other 
B — > PP decays, the theoretical predictions are still consistent with the measurements 
if one takes into account still large theoretical and experimental uncertainties. 

For penguin-dominated B — > PV decays, the new physics contributions in model III 
are significant, from 30% to 60%, and have a weak or moderate dependence on M H +, 
9 and N e " , as illustrated in the tables (|VTJJ-|VTrj) and figures (|n]-0). The CLEO 



measurements of BiB — > K* + rj, K* rj) are much larger than theoretical predictions in 
the SM and hence large new physics enhancements in model III are indeed helpful to 
lead to or improve the agreement between the data and theoretical predictions. 

8 In model III, the new physics contributions to different B — > VV decay modes are 
varying greatly: from —11% to ~ 110%. For decay modes B — > K* u, K* + (j), K*°cf), 
for example, the new physics enhancements are significant: ~ (60 — 110)%, and hence 
the theoretical predictions in model III are close to or slightly surpass the current 
CLEO upper limits. These decay modes will be observed soon. 
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APPENDIX A: INPUT PARAMETERS AND FORM FACTORS 



In this appendix we present relevant input parameters. We use the same set of input 
parameters for the quark masses, decay constants, Wolfenstein parameters and form factors 
as Ref. 0. 

• Input parameters of electroweak and strong coupling constant, gauge boson masses, 
B meson masses, light meson masses, • • •, are as follows (all masses in unit of GeV ) 



a em = 1/128, a s {M z ) = 0.118, sin 2 ^ = 0.23, G F = 1.16639 x 10-\GeV)~ 2 , 
M z = 91.187, M w = 80.41, m B o = m B ± = 5.279, m^± = 0.140, 
m n o = 0.135, m v = 0.547, ny = 0.958, m p = 0.770, m w = 0.782, 
m = 1.019, m K ± = 0.494, m K o = 0.498, m K *± = 0.892, m K *o = 0.896, 
r(Bi) = IMps, r(B° d ) = 1.56p«, (Al) 



• For the elements of CKM matrix, we use Wolfenstein parametrization, and fix the 
parameters A, A, p to their central values, A = 0.81, A = 0.2205, p = 0.12 and varying 
r\ in the range of rj = 0.34 ± 0.08. 

• We first treat the internal quark masses in the loops in connection with the function 
G{m) as constituent masses, 

m b = A.88GeV, m c = 1.5GeV, m s = 0.5GeV, m u = m d = 0.2GeV. (A2) 



Secondly, we will use the current quark masses for rrii (i = u, d, s, c, b) which appear 
through the equation of motion when working out the hadronic matrix elements. For 
fi = 2.5GeV, one finds 

m b = 4.88GeV, m c = l.BGeV, m s = 0.122GeV, m d = 7.6MeV, m u = A.2MeV. (A3) 



For the mass of heavy top quark we also use m t = mt(m t ) = lQ8GeV. 

For the decay constants of light mesons, the following values will be used in the nu- 
merical calculations (in the units of MeV): 

U = 133, f K = 158, f K . = 214, f p = 210, f u = 195, U = 233, 

/« = f d n = 78, % = = 68, f c n = -0.9, f^, = -0.23, 

f* = -113, ft = 141. (A4) 



where and / s (/) have been defined in the two- angle-mixing formalism with 6 = 
—9.1° and = —22.2° |Kj For more details about the mixings between r\ and rj', one 
can see [|5"0"| , |I5 |. 
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The form factors at the zero momentum transfer in the BSW model [1131 have been 



collected in Table 2 of Ref. 0. For the convenience of the reader we list them here: 

F B ^ n (0) = 0.33, F^ K (0) = 0.38, F B_n? (0) = 0.145, F B ^ r ''(0) = 0.135, 

<J(0) = = °- 28 ' A*-* K *(0) = 0.32, A B ? K *(0) = 0.33, 

V B ^ P (0) = V B ^(0) = 0.33, V B ^ K *(0) = 0.37, (A5) 

In the LQQSR approach, the form factors at zero momentum transfer being used in 
our numerical calculations are, 

F B_W (0) = 0.36, F B ^ K (0) = 0.41, F B ^ V (0) = 0.16, F B ^ r >' (0) = 0.145, 
{A , Ai, A 2 , V}(5 -»• p) = {0.30, 0.27, 0.26, 0.35}, 
{A , A u A 2 , V}(B -> iT) = {0.39, 0.35, 0.34, 0.48}, 

{A , Ai, A 2 , V}(B ^uj) = {0.30, 0.27, 0.26, 0.35}. (A6) 

The form factors F 0i i(fc 2 ), A 0t i^(k 2 ) and \^(fc 2 ) were defined in Ref. |Tj| as 

*b(0) F x (0) 



l-fc 2 /m 2 (0+)' iV ' l-k 2 /m 2 (l-)' 
M0) A A 1 (0) 



l-P/m 2 (0-)' v 7 l-P/m 2 (l+)' 

= y( * ) = i-»fti-) " (A7) 

The pole masses being used to evaluate the fc 2 -dependence of form factors are 

{m(0"), m(l"), m(l + ), m(0 + )} = {5.2789, 5.3248, 5.37, 5.73}, (A8) 

for ub and db currents, and 

{m(0"), m(l~), m(l + ), m(0 + )} = {5.3693, 5.41, 5.82, 5.89}, (A9) 

for sb currents. 
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TABLES 



TABLE I. Wilson coefficients d(fi) and C e J f {ji) in the SM and models I, II and III at the 
scale n = 2.5 GeV, with M H + = 200 GeV, tan/3 = 2 and 9 = 0°, 30°. 





SM 


Model I 


Model II 


Model III: 9 = 0° 


Model III: 9 = 30 


Ci 


1.1245 


1.1245 


1.1245 


1.1245 


1.1245 


c 2 


-0.2662 


-0.2662 


-0.2662 


-0.2662 


-0.2662 


c 3 


0.0186 


0.0187 


0.0187 


0.0186 


0.0186 


c A 


-0.0458 


-0.0458 


-0.0458 


-0.0458 


-0.0458 


c 5 


0.0113 


0.0113 


0.0113 


0.0113 


0.0113 


c 6 


-0.0587 


-0.0585 


-0.0585 


-0.0587 


-0.0587 


c 7 


0.0006 


0.0006 


0.0006 


0.0006 


0.0006 


c 8 


0.0007 


0.0007 


0.0007 


0.0007 


0.0007 


c 9 


-0.0095 


-0.0099 


-0.0099 


-0.0096 


-0.0096 


Cio 


0.0026 


0.0027 


0.0027 


0.0026 


0.0026 


C^ff 


-0.1527 


-0.1321 


-0.2487 


0.3364 


0.2708 + 0.2448i 
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TABLE II. Numerical values of for the transitions b — > d \b — > d ]. The first, second and 
third entries for 03, ... , a±o refer to the values of in the SM and models II and III, respectively. 
All entries for 03, ... , a 10 should be multiplied with 1CT 4 . 





N e ff = 2 


N e ff = 3 


N e ff = 00 


ai 


0.995 [0.995] 




1 nci M oftil 
l.Ubl [l.UblJ 


1.192 [1.192] 




0.201 [0.201] 




n nn 1 ? rn nn^l 


-0.395 [-0.395] 


^3 


1 CI 

— lb 


— [— zo — 16%\ 




77 f77l 

77 [77] 


261 + 13*' [280 + 47*'] 




— 10 


— 7z [—19 — 16i\ 




77 [77] 


252 + 13*' [271 + 47*'] 




a n 

—40 


— 1% [—49 — l6%\ 




77 T771 

77 [77] 


310 + 13*' [329 + 47*] 


04 


— 4zo - 


- 36i [—470 — 117zJ 


— 4b7 


oct'F n ^ 1 or '1 

— 35* [—517 — 125*] 


-554 - 39* [-610 - 141*] 




—398 - 


- 33* [—445 — 117*J 


—440 


oc; f 4 n n 1 r»r;l 

— dot [—490 — llbi\ 


-524 - 39*' [-581 - 141*] 




—54b - 


00,' r trno 1 1 ^7„"l 

- 33* [—592 — 117*j 


KOI 
— 0» 1 


— oDZ [— D4o — 1Z02J 


-701 - 39* [-757 - 141*] 


a 5 


1 oq 

— i»o 


7/." r ono OQ„'l 

— tl [— zUz — zo*J 




-71 [-71] 


171 + 13*' [190 + 47*] 




1 Q7 

— 187 


— (I [—190 — 24*J 




-71 [-71] 


161 + 13*' [180 + 47*] 




Ol 7 

—21 / 


7' r 00c 00 'l 
— 11 [— zzb — zo*J 




71 r 71 1 
-a [-ti-\ 


220 + 13*' [239 + 47*] 


a 6 


— b4z - 


00; f con 11 7i'l 

- oo* [—boy — lli*J 


£71 

—01 1 


or* f 701 1 OCC'l 

— 35* {—III — lzo*J 


-728 - 39*' [-784 - 141*] 




—bib - 


- oo* [— bbo — 11 ii\ 


ft A O 

— b4z 


— 35* [— by3 — 125*J 


-696 - 39*' [-752 - 141*] 




— /b4 - 


00 • r 01 1 11 7~1 
- 661 [—oil — 11 ll\ 


— oUl 


oka r oci 1 oc;„i 
— 00* [—oOl — 1ZDIJ 


-874 - 39* [-931 - 141*] 


a 7 


8.1 - 


-0.9? [7.7-1.7*] 


ft Q 

0.0 


n n„' rc a 1 7„'l 
— U.y* [b.4 — 1. 1%\ 


4.3 - 0.9* [3.9 - 1.7*] 




9.3 - 


-0.9* [8.9- 1.7*] 


q n 
o.U 


n OA T7 c 1 7„'l 

— u.y* [1.0 — 1. /*j 


5.3 - 0.9*' [4.9 - 1.7*] 




8.3- 


- 0.9*' [7.9-1.7*] 


7.0 


- 0.9* [6.6 - 1.7*] 


4.5 - 0.9* [4.1 - 1.7*'] 


a 8 


9.7- 


- 0.5* [9.5 - 0.8*] 


9.0 


- 0.3* [8.8 - 0.6*] 


7.5 [7.5] 




11 - 


- 0.5*' [11 - 0.8*] 


9.9 


- 0.3*' [9.7 - 0.6*] 


8.1 [8.1] 




9.9- 


- 0.5* [9.7 - 0.8*] 


9.1 


- 0.3* [9.0 - 0.6*] 


7.6 [7.6] 


a 9 


-84- 


-0.9* [-84-1.7*'] 


-90 


- 0.9*' [-90 - 1.7*] 


-102 - 0.9*' [-102 - 1.7*'] 




-87- 


- 0.9* [-87-1.7*] 


-93 


- 0.9*' [-94 - 1.7*] 


-106 - 0.9*' [-106 - 1.7*] 




-84- 


- 0.9* [-85 - 1.7*] 


-90 


- 0.9*' [-91 - 1.7*] 


-103 - 0.9*' [-103 - 1.7*] 




-14- 


-0.5* [-15-0.8*] 


2.6 


- 0.3* [2.5 - 0.6*] 


37 [37] 




-15- 


- 0.5* [-15 - 0.8*] 


2.8 


- 0.3*' [2.7 - 0.6*] 


38 [38] 




-15- 


-0.5* [-15-0.8*] 


2.7 


- 0.3* [2.5 - 0.6*] 


37 [37] 
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TABLE III. Same as Table but for b -> s [b->s] transitions. 





N e ff = 2 


N e ff = 3 




= 00 


Cl^ 




0.995 


0.995 






1.061 


[1.061] 


1.192 


"1.192] 


a>2 




0.201 


0.201 






0.026 


[0.026] 


-0.395 


-0.395] 




-21 


- 14? 


-19 - 


- 14?] 




77 


[77] 


272 + 29? 


"269 + 29?'] 




-15 


— 14? 


— 14 - 


- 14?J 




77 


[77] 


262 + 29? 


260 + 29?'] 




-45 


- 14? 


-44 - 


- 14?'] 




77 


[77] 


320 + 29? 


318 + 29?'] 


(24 


-449 


- 72? 


-442 


- 72?] 


-494 


- 77? 


[-487 - 77?'] 


-585 - 86? 


"_576 _ 86?'] 




—424 


— 72i 


—417 


- 72?'] 


-468 


- 77? 


[-460 - 77?] 


-555 - 87? 


[_547 _ 87?'] 




-571 


— 72i 


-564 


- 72?] 


-625 


- 77? 


[-617 - 77?] 


-732 - 86? 


[-723 - 86?] 


^0 


-198 


- 14? 


-196 


- 14?] 




-71 


[-71] 


181 + 29? 


"179 + 29?'] 




-192 


- 14? 


-191 


-1 A '1 

- 14?] 




-71 


r r—r 1 1 

[-71] 


172 + 29? 


[169 + 29?'] 




-222 


- 14? 


-221 


- 14?'] 




-71 


[-71] 


230 + 29? 


[228 + 29?'] 




-667 


- 72i 


-660 


- 72?] 


-698 


- 77? 


[-691 - 77?'] 


-758 - 86? 


"-750 - 86?'] 




-641 


- 72i 


-635 


- 72?'] 


-670 


- 77? 


[-663 - 77?'] 


-727 - 87? 


_719 _ 87i] 




-790 


- 72i 


-783 


- 72?'] 


-828 


- 77? 


[-821 - 77?'] 


-905 - 86? 


[-897 - 87?] 


(27 


7.9 - 


- 1.3? 


7.9 - 


1.3?] 


6.6 - 


- 1.3? 


[6.7 - 1.3?'] 


4.1 - 1.3? 


r 4.2 - 1.3?'] 




9.1 - 


- 1.3? 


9.2 - 


1.3?] 


7.7 - 


- 1.3? 


[7.8 - 1.3?'] 


5.0 - 1.3? 


7 5.1 - 1.3?'] 




8.1 - 


- 1.3? 


"8.2 - 


1.3?] 


6.8 - 


- 1.3? 


[6.9 - 1.3?'] 


4.3 - 1.3? 


[4.3 - 1.3?] 


(2s 


9.6 - 


- 0.6? 


9.6 - 


0.6?'] 


8.9 - 


- 0.4? 


[8.9 - 0.4?'] 


7.5 


[7.51 

L J 




10.6 - 


- 0.6? 


10.6 - 


- 0.6?'] 


9.8 - 


- 0.4? 


[9.8 - 0.4?] 


8.1 


[8 1] 




9.8- 


- 0.6? 


9.8 - 


0.6?'] 


9.1 - 


-0.4?' 


[9.1 - 0.4?'] 


7.6 


[7.6] 


a 9 


-84- 


- 1.3? 


-84- 


- 1.3?'] 


-90 - 


- 1.3? 


[-90 - 1.3?] 


-102 - 1.3? 


[-102 - 1.3?'] 




-87- 


- 1.3?' 


-87- 


- 1.3?] 


-94- 


- 1.3? 


[-94-1.3?] 


-106 - 1.3?' 


-106 - 1.3?'] 




-85 - 


- 1.3* 


-84- 


- 1.3?'] 


-91 - 


- 1.3? 


[-91 - 1.3?] 


-103 - 1.3? 


[-103 - 1.3?'] 




-15 - 


- 0.6?' 


-14- 


- 0.6?] 


2.6 - 


-0.4? 


[2.6 - 0.4?] 


37 


[37] 




-15 - 


- 0.6?' 


-15 - 


- 0.6?'] 


2.8 - 


- 0.4? 


[2.8 - 0.4?] 


38 


[38] 




-15 - 


- 0.6?' 


-15 - 


- 0.6?] 


2.6 - 


-0.4i 


[2.6 - 0.4?] 


37 


[37] 
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TABLE IV. B{B -► PP) (in units of 10~ 6 ) in the SM using the BSW [LQSSR] form factors, 
with k 2 = ml/2 and N e ^ = 2, 3, oo. The last column shows the CLEO measurements and upper 
limits at 90%C.L. [22-25]. 



Channel 


Class 


N e ff = 2 


N e ff = 3 


N e ff = oo 


Data 


B° -► 7T+7T- 

B° - 7r°7r° 
- vr+vr 


I 

II 
III 


9.10 [10.8] 
0.28 [0.33] 
6.41 [7.62] 


10.3 [12.3] 
0.15 [0.18] 
5.06 [6.02] 


13.0 [15.5] 
0.92 [1.09] 
2.85 [3.39] 


4.3i^±0.5 

< 9.3 

< 12.7 


5° — ► 7/7/ 
B° -► 7/7/' 
5° 7/' 7]' 


II 
II 
II 


0.14 [0.17] 
0.14 [0.17] 
0.04 [0.05] 


0.10 [0.13] 
0.08 [0.09] 
0.01 [0.01] 


0.29 [0.36] 
0.38 [0.45] 
0.13 [0.15] 


< 18 

< 27 

< 47 


B+ _> vr+r; 
5+ -► 7T+7/' 
5° - 7T°7/ 
5° TT^Tj' 


III 
III 

V 
V 


3.51 [4.25] 
2.49 [2.90 
0.26 [0.31] 
0.06 [0.07] 


2.78 [3.37] 
1.88 [2.17] 
0.29 [0.35] 
0.08 [0.09] 


1.75 [2.13] 
1.02 [1.17] 
0.39 [0.47] 
0.14 [0.17] 


< 5.7 

< 12 

< 2.9 

< 5.7 


B+ 

5 + tf 7T+ 
£° - K 7T° 


IV 
IV 
IV 
IV 


12.0 [14.3] 

17.8 [21.2] 

19.9 [23.7] 
7.27 [8.68] 


13.5 [16.0] 
19.8 [23.5] 
23.2 [27.7] 
8.31 [9.92] 


16.7 [19.8] 
24.0 [28.5] 

30.6 [36.4] 

10.7 [12.7] 


n 5+3.0+1.4 

17.2±|;| ± 1.2 
18.2±|-g ± 1.6 

14 6+5^ +2.4 
14 - D -5.1 -3.3 


5+ -► K+7? 
5+ -»• K+7/' 
5° -► K°7/ 
5° K°7/' 


IV 
IV 
IV 
IV 


3.91 [4.37] 
22.6 [26.2] 
3.22 [3.57] 
21.9 [25.5] 


4.56 [5.10] 
28.5 [33.1] 
3.63 [4.02] 
28.2 [32.7] 


6.07 [6.80] 
42.4 [49.2] 
4.58 [5.07] 
43.0 [49.9] 


< 6.9 
80+^° ± 7 

< 9.3 
89±i| ± 9 


B ° - ^0^0 


IV 
IV 


1.16 [1.35] 
1.10 [1.28] 


1.35 [1.58] 
1.28 [1.49] 


1.78 [2.07] 
1.68 [1.96] 


< 5.1 

< 17 
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TABLE V. B(B -> PP) (in units of 1(T 6 ) in model III using the BSW form factors, with 
k 2 = ml/2, N e M = 2, 3, oo, M H+ = 200GeV and 9 = 0°,30°, respectively. 





= 0° 


5B [%] 


= 30° 


5B [%] 


Channel 


2 


3 


00 


2 


3 


00 


2 


3 


00 


2 


3 


00 


B U — ► 7T + 7T 


9.33 


10.6 


13.3 


2.5 


2.5 


2.4 


8.83 


10.0 


12.6 


-3.0 


-3.1 


-3.1 


dO 
i? u — > 7T U 7T U 


0.36 


0.25 


1.03 


30 


61 


13 


0.39 


0.23 


0.92 


40 


52 


—0.5 


— > 7r^7r 


6.41 


5.06 


2.85 


0.0 


0.0 


0.0 


6.41 


5.06 


2.85 


0.0 


0.0 


0.0 


B — > 7777 


0.18 


0.15 


0.36 


29 


47 


21 


0.16 


0.14 


0.38 


15 


39 


29 


-B U — > T/T/ 


0.19 


0.13 


0.46 


29 


68 


20 


0.16 


0.12 


0.506 


9.6 


57 


30 


5 — > 77 77 


0.05 


0.02 


0.15 


20 


127 


17 


0.04 


0.02 


0.16 


—5.4 


103 


30 


D+ + 


3.82 


3.13 


2.20 


8.7 


13 


26 


3.48 


2.81 


1.95 


-1.1 


1.2 


12 


7^)4- -1- / 

i? + — ► 7r + r/ 


2.63 


2.05 


1.27 


5.4 


9.0 


24 


2.37 


1.81 


1.09 


-4.8 


-3.5 


6.5 


nn o 
B u — > 7r u r/ 


0.39 


0.44 


0.59 


50 


51 


49 


0.36 


0.42 


0.57 


39 


43 


46 


-B U — > 7T U 77 


0.11 


0.14 


0.25 


92 


91 


72 


0.10 


0.13 


0.25 


65 


76 


70 


D+ 7^ + 


17.4 


19.6 


24.4 


45 


45 


46 


17.2 


19.3 


23.8 


43 


43 


43 


7^0 7^4- — 


26.8 


29.9 


36.5 


51 


51 


53 


26.5 


29.5 


36.1 


49 


49 


51 


5 + -► K°7T+ 


29.8 


34.6 


45.3 


50 


49 


48 


28.6 


33.3 


43.5 


44 


43 


42 


£° - K°7T 


11.4 


13.0 


16.7 


57 


57 


56 


10.8 


12.5 


16.1 


49 


50 


51 


B+ K+r/ 


5.69 


6.63 


8.78 


45 


45 


45 


5.30 


6.22 


8.34 


36 


36 


37 


73+ _^ x+r/ 


38.0 


46.9 


67.5 


68 


65 


59 


36.8 


45.2 


64.9 


63 


59 


53 


5° -► iY°r/ 


4.86 


5.50 


6.93 


51 


51 


51 


4.63 


5.27 


6.73 


44 


45 


47 


5° iY°r/ 


36.7 


45.9 


67.3 


67 


63 


57 


35.1 


43.8 


64.2 


60 


56 


49 




1.73 


2.01 


2.62 


49 


48 


47 


1.64 


1.91 


2.49 


41 


41 


40 




1.64 


1.90 


2.48 


49 


48 


47 


1.55 


1.80 


2.36 


41 


41 


40 
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TABLE VI. B(B -> PP) (in units of 10" 6 ) in models I and II using the BSW form factors, 
with k 2 = ml/2, N e f-f = 2, 3, oo, tan/3 = 2 and M H + = 200GeV. 





Model I 


SB [%] 


Model II 


5B [%] 


Channel 


2 


3 


00 


2 


3 


00 


2 


3 


00 


2 


3 


00 


73° -»■ 7T+7T~ 


9.11 


10.3 


13.0 


0.1 


0.1 


0.1 


9.1 


10.3 


13.0 


-0.5 


-0.4 


-0.4 


5° -> 7r 7T° 


0.28 


0.15 


0.92 


—0.1 


—0.1 


0.1 


0.3 


0.1 


0.9 


—6.2 


—12.6 


—2.5 


n4- -1- 
_g+ — > 7T + 7T U 


6.41 


5.06 


2.85 


0.0 


0.0 


0.0 


6.4 


5.1 


2.9 


0.0 


0.0 


0.0 


/3 — > 7777 


0.14 


0.11 


0.30 


1.5 


2.4 


1.1 


0.1 


0.1 


0.3 


-4.7 


-7.5 


-3.4 


73° -> 7777' 


0.14 


0.08 


0.38 


1.0 


2.3 


0.7 


0.1 


0.1 


0.4 


-4.8 


-11.4 


-3.5 


B° -> 7/7/ 


0.04 


0.01 


0.13 


0.3 


2.0 


0.2 


0.04 


0.01 


0.1 


—2.7 


—20.2 


—3.1 


73+ — »■ 7T+7? 


3.53 


2.79 


1.77 


0.4 


0.5 


1.0 


3.5 


2.7 


1.7 


-1.5 


-2.1 


-4.40 


i3+ -> vr+r/ 


2.50 


1.88 


1.03 


0.1 


0.2 


0.4 


2.5 


1.9 


1.0 


-0.9 


-1.6 


-4.5 


B° -> 7r°r/ 


0.26 


0.30 


0.40 


1.5 


1.6 


1.7 


0.2 


0.3 


0.4 


-9.0 


-9.1 


-8.6 


i3° -> 7T 77 


0.06 


0.08 


0.15 


0.9 


0.9 


0.7 


0.05 


0.1 


0.1 


—15.9 


— 16.4 


— 13.5 


i3 + -> K+7T 


12.3 


13.8 


17.1 


2.4 


2.3 


2.2 


11.2 


12.5 


15.4 


-7.1 


-7.3 


-7.5 


73° -»■ 


18.0 


20.0 


24.2 


1.3 


1.3 


1.2 


16.2 


17.9 


21.6 


-9.3 


-9.4 


-9.7 


73 + K°7T+ 


20.2 


23.6 


31.1 


1.4 


1.5 


1.6 


18.1 


21.2 


28.0 


-9.0 


-8.8 


-8.5 


73° - K°7T 


7.28 


8.33 


10.7 


0.1 


0.2 


0.5 


6.4 


7.4 


9.5 


-11.6 


-11.5 


-11.1 


73+ -► K+r? 


3.91 


4.58 


6.12 


0.0 


0.3 


0.8 


3.5 


4.1 


5.6 


-9.6 


-9.3 


-8.7 


5+ -► K+r]' 


23.0 


28.9 


43.0 


1.6 


1.5 


1.4 


19.8 


25.2 


37.8 


-12.1 


-11.6 


-10.8 


B° -► i^°7? 


3.24 


3.66 


4.63 


0.6 


0.8 


1.1 


2.9 


3.3 


4.1 


-9.9 


-9.8 


-9.5 




22.3 


28.5 


43.5 


1.4 


1.4 


1.3 


19.3 


24.9 


38.5 


-12.2 


-11.5 


-10.5 


73+ K+i?° 


1.18 


1.37 


1.81 


1.4 


1.5 


1.6 


1.1 


1.2 


1.6 


-8.9 


-8.7 


-8.4 




1.11 


1.30 


1.71 


1.4 


1.5 


1.6 


1.0 


1.2 


1.5 


-8.9 


-8.7 


-8.4 
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TABLE VII. B -> PV branching ratios (in units of 10" 6 ) using the BSW [ LQQSR ] form 
factors in the SM, with with k 2 = ml/2, N e ff = 2, 3, oo. The last column shows the CLEO 
measurements and upper limits (90% C.L.) [22-25,27]. 



Channel 


Class 


Are// = 2 


N e ff = 3 


N e Jf = oo 


Data 


B" - p+TT- 
B° P -TT+ 

B° - p°7r° 

B+ p°7T+ 

B + p+vrO 


I 
I 

II 

III 

III 


21.1 [25. 1J 
5.7 [6.5] 
0.49 [0.58] 
5.72 [6.63] 
13.5 [16.0] 


24.0 [28.5] 

6.5 [7.4] 
0.06 [0.07] 
3.46 [3.97] 
12.6 [15.0] 


30.3 [36.0] 

8.2 [9.4] 
2.05 [2.41] 
0.71 [0.78] 
10.9 [13.1] 


1 97 fi+ 8 ' 4 -1-4 9 

< 5.1 

10.4l^ ± 2.1 

< 43 


5° -► p°r7 
5° pV 
5+ -> p+r/ 
B+ p+r/ 


II 
II 
III 
III 


0.01 [0.02] 
0.01 [0.01] 
5.44 [6.57] 
4.35 [5.02] 


0.02 [0.02] 
0.002 [0.003] 
4.75 [5.79] 
3.81 [4.40] 


0.06 [0.08] 
0.03 [0.03] 
3.54 [4.38] 
2.85 [3.29] 


< 10 

< 12 

< 15 

< 33 


5+ -► o;vr+ 
B -► wry 
5° -► wr/ 


II 
III 
II 
II 


0.29 [0.35] 
6.32 [7.35] 
0.32 [0.38] 
0.20 [0.23] 


0.08 [0.09] 
3.75 [4.31] 
0.03 [0.04] 
0.001 [0.002] 


0.15 [0.19] 
0.78 [0.85] 
0.82 [0.98] 
0.68 [0.79] 


< 5.5 

11.3j:l;| ± 1A 

< 12 

< 60 


B° 0vr° 
5° -► 07/ 


V 
V 
V 
V 


0.03 [0.04] 
0.06 [0.08] 
0.01 [0.01] 
0.01 [0.01] 


0.002 [0.002] 
0.004 [0.005] 
0.001 [0.001] 
0.001 [0.001] 


0.23 [0.27] 
0.49 [0.58] 
0.09 [0.10] 
0.07 [0.08] 


< 5.4 

< 4 

< 9 

< 31 


B o _^ k*°K° 
B + _> 


IV 
IV 
V 
IV 


0.42 [0.49] 
0.40 [0.46] 
0.005 [0.007] 
0.004 [0.006] 


0.53 [0.61] 
0.50 [0.58] 
0.002 [0.003] 
0.002 [0.003] 


0.78 [0.90] 
0.73 [0.89] 
0.001 [0.001] 
0.001 [0.001] 


< 5.3 

< 12 


5° p°K° 

B + ^ p0 K + 

B° -► p-K+ 
B + ^ p+^0 


IV 
IV 

I 

IV 


0.52 [0.60] 
0.39 [0.46] 
0.54 [0.62] 
0.11 [0.12] 


0.53 [0.62] 
0.31 [0.36] 
0.59 [0.68] 
0.05 [0.05] 


0.71 [0.83] 
0.31 [0.36] 
0.70 [0.81] 
0.01 [0.01] 


< 27 

< 17 

< 25 

< 48 


5+ -► K*+r/ 
5+ -► K*+r/ 
5° -► K* r? 
B° K*%' 


IV 

III 

IV 

V 


2.43 [3.12] 
0.66 [1.14] 
2.32 [2.98] 
0.33 [0.69] 


2.39 [3.04] 
0.36 [0.61] 
2.54 [3.23] 
0.09 [0.23] 


2.32 [2.89] 
0.24 [0.23] 
3.06 [3.82] 
0.31 [0.26] 


26.4l^ ± 3.3 

< 35 

13.8l!;jj ± 1-6 

< 20 


Y J) — r i V 71 

B o ^ K * 0vr o 


TV 
1 V 

IV 
IV 
IV 


r ri n 9i 
o.oy [iu.zj 

2.44 [2.77] 

4.95 [6.09] 

7.35 [8.75] 


3.02 [3.43] 
5.55 [6.84] 
9.23 [11.0] 


1Z.U [11. Ol 

4.42 [5.01] 
6.91 [8.52] 
13.6 [16.2] 


99+8 +4 
zz -6-5 

< 3.6 

< 31 

< 16 


5+ -► 
5° -► 


V 
V 


22.1 [25.7] 
20.9 [24.3] 


11.5 [13.4] 
10.9 [12.6] 


0.60 [0.70] 
0.57 [0.66] 


17.2+|;^ ± 1.8 
< 28 


5° -► 

5+ -► wK+ 


V 
V 


3.31 [3.86] 
3.53 [4.11] 


0.002 [0.003] 
0.25 [0.28] 


13.3 [15.4] 
16.5 [19.2] 


< 21 

< 7.9 
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TABLE VIII. B -> PV branching ratios (in units of 1(T 6 ) using the BSW form factors in 
model III, assuming M H + = 200GeV, 9 = 0° and N e ^ = 2, 3, oo. 







SM 


Model III 


5B [%} 


Channel 


Class 


2 


3 


oo 


2 


3 


oo 


2 


3 


oo 


nf) -4- 

B" — »■ /9+7T 


I 


21.1 


24.0 


on o 

30.3 


O "I O 

21.2 


24.1 


on r* 

30.5 


n t 

0.7 


n t 

0.7 


n >"7 

0.7 


B v — ► p 7T + 


I 


5.70 


6.48 


8.19 


5.70 


6.48 


8.19 


0.0 


0.0 


0.0 


r>o n n 

B U — ► /9 U 7T U 


II 


0.49 


0.06 


2.05 


0.54 


0.11 


2.12 


9.8 


99.6 


3.5 


T-i-t- n 4- 

5 + — »■ p U 7T + 


III 


5.72 


3.46 


0.71 


5.79 


3.54 


0.81 


1.3 


2.3 


14.0 


7->^ -1- 

i> + — > p + 7T U 


III 


1 9 K 

13.5 


12.0 


inn 
10.9 


13.6 


1 O T 

12. 7 


i 1 n 
11.0 


n A 
0.4 


n k 
0.5 


n ^7 
0. / 


D n 


TT 

11 


0.01 


0.02 


0.06 


0.03 


0.03 


0.08 


86.0 


100 


40.1 


r^n n / 
£r — > /9 U ?7 


II 


0.01 


0.003 


0.03 


0.004 


0.001 


0.03 


—47.3 


—54.4 


18.1 


n-L -U 

5+ — > p+ry 


III 


5.44 


4.75 


3.54 


5.46 


4.79 


3.59 


0.5 


0.7 


1.4 


£? + — > p+57' 


III 


A 9 £ 


3.81 


2.85 


A OA 

4.34 


9 01 

3.81 


2.86 


n o 
— U.2 


n no 
— U.U8 


n A 
0.4 


dO 


TT 

11 


0.29 


0.08 


0.15 


0.45 


0.14 


0.15 


54.4 


77.0 


0.8 


n-l- -i- 
B + — >■ L07T 


III 


6.32 


3.75 


0.78 


6.63 


3.86 


0.79 


5.0 


3.1 


1.1 


B" — > cj?7 


II 


0.32 


0.03 


0.82 


0.37 


0.05 


0.83 


16.3 


67.1 


0.2 


7^>n / 
B u — > wr/' 


II 


0.20 


n nm 

0.001 


0.68 


n oo 

0.22 


0.004 


0.69 


n c 

9.5 


155 


n cr 

0.5 


n0 i 


V 


0.03 


0.002 


0.23 


0.05 


0.002 


0.33 


59.1 


1.9 


42.3 


7">4- i 4- 
±> + — > 07T + 


V 


0.06 


0.004 


0.49 


0.10 


0.004 


0.69 


59.1 


1.9 


42.3 


£? U — ► 077 


V 


0.01 


0.001 


0.09 


0.02 


0.001 


0.12 


59.1 


1.9 


42.3 


nf) / / 

£? U — ► (/)7/ / 


V 


n m 
0.01 


n nm 
0.001 


n n^7 
0.0 7 


n m 
0.01 


n nm 
0.001 


n i n 
0.10 


59.1 


i n 
1.9 


O 9 

42.3 


n-l- TV*0 TV-l- 

B^ — > K J K^ 


TT T 

IV 


0.42 


0.53 


0.78 


0.68 


0.83 


1.19 


61.0 


57.9 


53.3 


B v — > K K 


IV 


0.40 


0.50 


0.73 


0.64 


0.79 


1.12 


61.0 


57.9 


53.3 


T"~» — I— TV*— 1— 7" t 


V 


0.005 


0.002 


0.001 


0.002 


0.001 


0.003 


-58.4 


-72.5 


256 


— > if* u if u 


IV 


n nn a 

0.004 


n nno 

0.002 


n nn i 

0.001 


0.002 


0.001 


n nno 

0.003 


—58.4 


— 72.5 


256 


n0 jy-0 
±r — > p u A u 


IV 


0.52 


0.53 


0.71 


0.43 


0.44 


0.60 


-17.1 


-18.3 


-16.5 


n-l- tv-I- 

5+ — »■ p u if + 


IV 


0.39 


0.31 


0.31 


0.43 


0.36 


0.40 


8.0 


16.4 


30.6 


nD — TV- 4- 

£r — > p if + 


I 


0.54 


0.59 


0.70 


0.47 


0.52 


0.62 


-13.1 


— 12.7 


— 11.8 


n-l- -I- Tvfl 

5+ — > p+if u 


IV 


nil 
0.11 


0.05 


n m 
0.01 


n n £ 
0.05 


n m 
0.01 


n no 
0.02 


pen o 
—50.8 


^7n i 
— /0.1 


9£?9 

363 


J ") -1- TV" * -1- 

— > K T] 


TT T 

IV 


2.43 


2.39 


2.32 


3.27 


3.29 


3.34 


34.4 


37.4 


43.5 


n-l- TV*-!— / 

— ► K ^Tj 


III 


0.66 


n or 
0.36 


n i A 
U.24 


n 9 1 
0.31 




0.65 


ECO O 

— 02.2 


9 A O 

— 34.2 




n0 7y * ^ 

±> — ► A 77 


IV 


2.32 


2.54 


3.06 


3.15 


3.47 


4.20 


35.8 


36.5 


OT /( 

37.4 


nO 7V-*fl / 

i? u — > K 7] 


V 


0.33 


0.09 


0.31 


0.08 


0.10 


0.96 


-77.3 


6.9 


204 




IV 


8.59 


9.67 


12.0 


13.6 


15.4 


19.1 


58.6 


58.8 


59.3 


5 o ^ K * 0vr o 


IV 


2.44 


3.02 


4.42 


4.26 


5.18 


7.34 


74.9 


71.6 


66.0 


5+ if * + 7T° 


IV 


4.95 


5.55 


6.91 


7.42 


8.38 


10.5 


49.9 


50.9 


52.2 




IV 


7.35 


9.23 


13.6 


11.9 


14.7 


21.0 


62.0 


58.8 


54.1 


5+ -► 0if + 


V 


22.1 


11.5 


0.60 


35.7 


19.0 


1.29 


61.5 


65.3 


113 


5° - 0if 


V 


20.9 


10.9 


0.57 


33.7 


18.0 


1.21 


61.5 


65.3 


113 


B° - a;if 


V 


3.31 


0.002 


13.3 


5.33 


0.01 


19.4 


60.9 


175 


46.5 


5+ -> Cjif + 


V 


3.53 


0.25 


16.5 


5.57 


0.23 


23.6 


57.8 


-7.7 


42.9 
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TABLE IX. B — > branching ratios (in units of 10 6 ) using the BSW form factors in models 
I and II, assuming M H + = 200GeV, tan/3 = 2 and iV e// = 2, 3, oo. 





Model I 


SB [% 




Model II 


SB [%} 


Channel 


2 


3 


oo 


2 


3 


oo 


2 


3 


oo 


2 


3 


oo 


B° p+TT- 


21.1 


o i n 

24.0 


on o 

30.3 


0.0 


0.0 


0.0 


21.1 


oo n 

23.9 


on o 

30.2 


n i 

—0.1 


n i 

—0.1 


n i 

—0.1 


£° -► /9-7T+ 


5.70 


6.48 


8.19 


—0.0 


-0.0 


0.0 


5.70 


6.48 


8.19 


0.0 


0.0 


0.0 


£° - A° 


0.49 


0.06 


2.05 


-0.1 


-0.4 


0.0 


0.48 


0.05 


2.04 


—1.9 


-18.9 


-0.6 




5.72 


3.46 


0.71 


-0.0 


-0.0 


0.1 


5.70 


3.44 


0.69 


—0.3 


—0.5 


—2.8 


5 + p+TT 


1 9 K 

13.5 


12.0 


inn 
10.9 


n n 
0.0 


n n 
0.0 


0.0 


1 9 £ 

13.5 


in /; 
12.0 


inn 
10.9 


n i 
—0.1 


n i 
—0.1 


n i 
—0.1 


5° -► p 7? 


0.02 


0.02 


0.06 


2.9 


3.7 


1.6 


0.01 


0.01 


0.05 


— 12.8 


— 15.2 


—6.2 


5° A' 


0.01 


0.002 


0.03 


-0.2 


-0.1 


0.0 


0.01 


0.003 


0.03 


18.1 


39.1 


—1.1 


5+ -»• p+r/ 


5.44 


4.76 


3.54 


0.0 


0.0 


0.1 


5.43 


4.75 


3.53 


—0.1 


—0.1 


—0.2 


5+ -► p+r/ 


A 9 £ 


3.81 


2.85 


n n 
—0.0 


n n 

— u.u 


0.0 


A OP 

4.36 


9 Q 1 

3.81 


2.85 


n 1 
0.1 


n 1 
0.1 


n n 
—0.0 


5° - W 7T° 


0.29 


0.08 


0.15 


0.2 


-0.4 


0.5 


0.26 


0.07 


0.15 


-10.7 


— 15.4 


0.4 


5 + -> + 


6.33 


3.75 


0.78 


0.1 


0.1 


0.0 


6.26 


3.73 


0.78 


—0.9 


—0.5 


—0.2 




0.32 


0.03 


0.82 


n *i 

0.7 


3.6 


0.0 


0.31 


0.03 


0.82 


O 'I 

—2.7 


—9.9 


—0.1 


£° -► urf 


n onn 

0.209 


0.001 


0.68 


C\ o 

0.2 


n o 

0.3 


0.1 


n on 

0.20 


n nno 

0.002 


n r'O 

0.68 


1 O 

— 1.3 


8.7 


n i 

—0.1 


B° - 07T° 


0.03 


0.002 


0.23 


-0.8 


10.9 


2.4 


0.03 


0.002 


0.21 


— 12.9 


10.9 


-6.7 


5+ -»• 07T+ 


0.06 


0.004 


0.50 


—0.8 


10.9 


2.4 


0.06 


0.005 


0.45 


— 12.9 


10.9 


-6.7 


5° -> 07? 


0.01 


0.001 


0.09 


—0.8 


10.9 


2.4 


0.01 


0.001 


0.08 


— 12.9 


10.9 


-6.7 


5° -► 0T?' 


n m 
0.01 


n nm 
0.001 


0.0 7 


n o 
—0.8 


inn 
10.9 


2.4 


n ni 
0.01 


n nm 
0.001 


n r\p 
0.06 


ion 

— 12. y 


inn 
10.9 


—6. 1 


5+ _^ K*°K+ 


0.43 


0.54 


0.79 


1.8 


1.9 


2.1 


0.37 


0.47 


0.71 


— 10.7 


— 10.0 


-9.0 


B o _^ k*°K° 


0.40 


0.51 


0.75 


1.8 


1.9 


2.1 


0.35 


0.45 


0.67 


-10.7 


-10.0 


-9.0 


B + k*+K° 


0.005 


0.002 


0.001 


-3.5 


-7.1 


13.5 


0.01 


0.002 


0.001 


14.3 


20.9 


-16.6 


B o _^ k*°K° 


r\ r\r\ A 

0.004 


n nno 

0.002 


n nn i 

0.001 


—3.5 


—7.1 


13.5 


n nn c 

0.005 


n nno 

0.002 


n nni 

0.001 


1 A O 

14.3 


on n 

20.9 


— 16.6 


B° - p°K° 


0.53 


0.55 


0.74 


3.1 


3.3 


3.1 


0.56 


0.58 


0.77 


7.6 


8.2 


7.6 


B + ^ p o K + 


0.40 


0.32 


0.34 


2.4 


4.9 


9.0 


0.40 


0.31 


0.32 


0.7 


1.4 


2.4 




0.53 


0.58 


0.70 


—1.8 


— 1.4 


-0.6 


0.55 


0.60 


0.72 


1.6 


2.0 


2.7 


5+ p+K° 


n i n 
0.10 


n ri/i 
0.04 


n nn k 
0.005 


— o.U 


P 9 

—0.3 


11.6 


n i o 
U.12 


n nc; 
0.05 


n nn^ 
0.005 


1 1 s 

11.0 


1 P. T 

16. 1 


—6.5 


5+ -► K* + 7/ 


2.49 


2.45 


2.36 


2.5 


2.3 


1.7 


2.31 


2.26 


2.15 


—4.9 


—5.8 


—7.5 


5 + K* + r/ 


0.64 


0.35 


n o A 
U.24 


O /I 

— Z.4 


1 9 

—2.3 


2.0 


n 

0.1 1 


n A 9 
0.43 


n o 1 
U.21 


Ib.b 


1 n 

19.5 


1 o cc 

— Iz.o 


5° -► K*°77 


2.37 


2.60 


3.13 


2.4 


2.3 


2.2 


2.19 


2.40 


2.88 


—5.4 


—5.6 


r n 

—5.9 


5° -► £T*V 


0.32 


0.09 


0.33 


-3.3 


-4.7 


5.0 


0.42 


0.13 


0.23 


27.3 


42.6 


-25.4 


5° -► K* + 7T- 


8.76 


9.84 


12.2 


2.0 


1.8 


1.4 


7.74 


8.69 


10.8 


-9.8 


-10.1 


-10.5 


#o ^ K * 0vr o 


2.44 


3.03 


4.45 


0.1 


0.4 


0.8 


2.08 


2.60 


3.87 


-14.7 


-13.9 


-12.5 




5.11 


5.72 


7.10 


3.2 


3.0 


2.7 


4.60 


5.14 


6.35 


-7.2 


-7.5 


-8.1 




7.48 


9.41 


13.9 


1.8 


1.9 


2.1 


6.55 


8.29 


12.4 


-10.8 


-10.2 


-9.2 


5+ -► 


22.3 


11.6 


0.61 


1.1 


1.0 


0.7 


19.5 


10.1 


0.48 


-11.5 


-12.2 


-19.9 


B° - 


21.1 


11.0 


0.57 


1.1 


1.0 


0.7 


18.5 


9.54 


0.46 


-11.5 


-12.2 


-19.9 


5° - 


3.35 


0.002 


13.4 


1.2 


1.1 


1.3 


2.94 


0.003 


12.1 


-11.3 


25.7 


-8.5 


5+ -► 


3.59 


0.25 


16.7 


1.6 


-0.1 


1.4 


3.17 


0.25 


15.2 


-10.1 


2.3 


-7.7 
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TABLE X. B -> VV branching ratios (in units of 1CT 6 ) using the BSW [ LQQSR ] form factors 
in the SM, with k 2 = ml/2, N e ^ = 2, 3, oo. The last column shows the CLEO upper limits (90% 
C.L.) [22-25]. 



Channel 


Class 


N e ff = 2 


Are// = 3 


fyfeff 


= oo 


Data 


B° p+p- 


I 


17.8 


[19.8] 


20.2 [22.5] 


25.5 


[28.4] 


< 2200 


B° -► p°p° 


II 


0.39 


[0.44] 


0.09 [0.10] 


1.56 


[1.73] 


< 4.8 


B° _> coco 


II 


0.81 


[0.90] 


0.15 [0.17] 


1.22 


"1.35] 


< 19 


B+ 


III 


12.8 


[14.3] 


10.1 [11.3] 


5.69 


[6.33] 


< 120 




III 


15.7 


[17.4] 


12.2 [13.5] 


6.69 


[7.45] 


< 47 


#o _^ K* + p~ 


IV 


6.17 


[6.82] 


6.95 [7.68] 


8.64 


[9.55] 






IV 


1.73 


[1.82] 


2.01 [2.11] 


2.79 


[2.91] 


< 16.1 




IV 


5.22 


[5.97] 


5.97 [6.82] 


7.76 


[8.91] 


< 52 




IV 


6.65 


[7.35] 


8.36 [9.24] 


12.4 


[13.7] 






IV 


0.38 


[0.49] 


0.48 [0.61] 


0.70 


0.90] 


< 62 


B° — ► i^* J^*° 


IV 


0.36 


[0.47] 


0.46 [0.58] 


0.67 


[0.86] 


< 7.4 


5° - p°u; 


V 


0.45 


[0.50] 


0.24 [0.27] 


0.02 


[0.02] 


< 11 


5° -► K*°u; 


V 


13.5 


[16.1] 


4.52 [5.03] 


1.04 


[1.79] 


< 19 




V 


13.4 


[16.1] 


3.94 [4.39] 


2.74 


[4.01] 


< 52 


£ + -► + 


V 


21.8 


[27.8] 


11.3 [14.5] 


0.60 


[0.76] 


< 41 


5° -► k*V 


V 


20.6 


[26.2] 


10.7 [13.6] 


0.56 


[0.72] 


< 21 


5+ -► p+0 


V 


0.06 


[0.07] 


0.004 [0.005] 


0.47 


[0.52] 


< 16 


B° -► 


V 


0.03 


[0.03] 


0.001 [0.002] 


0.22 


[0.25] 


< 13 


B° ^ujcj) 


V 


0.03 


[0.03] 


0.001 [0.002] 


0.22 


[0.24] 


< 21 
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TABLE XI. B — ► VV branching ratios (in units of 10 6 ) using the BSW form factors in model 
III, assuming M H + = 200 GeV, N e ^ = 2, 3, oo and 9 = 0°, 30°. 



Channel 


9 = 0° 


5B [%} 


9 = 30° 


SB [%} 




2 


3 


oo 


2 


3 


oo 


2 


3 


oo 


2 


3 


oo 


n(l -1- — 
B U — > /9 + /9 


17.9 


20.3 


25.7 


0.7 


0.7 


0.7 


19.3 


19.7 


24.8 


-2.6 


-2.6 


-2.6 


nfi n n 


0.46 


0.16 


1.65 


16 


79 


5.7 


0.56 


0.15 


1.50 


28 


67 


-4.0 


£? U — ► UJUJ 


1.05 


0.24 


1.23 


29 


54 


1.0 


1.02 


0.22 


1.19 


14 


45 


-2.3 




12.8 


10.1 


5.69 


0.0 


0.0 


0.0 


14.3 


10.1 


5.69 


0.0 


0.0 


0.0 


-4- 

— ► p+u; 


13.8 


10.9 


6.12 


-12 


-11 


-8.5 


15.2 


10.9 


6.12 


-13 


-11 


-8.5 


£? U — > if * + /9 


9.78 


11.0 


13.8 


59 


59 


59 


11.1 


11.3 


14.1 


62 


63 


63 


£r — > K "p u 


3.18 


3.70 


4.98 


84 


84 


78 


3.14 


3.50 


4.90 


72 


74 


75 


23 + — > K + p 


7.61 


8.75 


11.4 


46 


47 


47 


8.79 


8.80 


11.3 


47 


48 


45 


— > K u p + 


10.8 


13.3 


19.0 


62 


59 


54 


11.4 


12.7 


18.2 


55 


52 


48 


~r~i —I— Ty*-I— 

B + — > K*+_fC* u 


0.61 


0.75 


1.08 


61 


58 


53 


0.73 


0.71 


1.02 


51 


49 


45 




0.59 


0.72 


1.03 


61 


58 


53 


0.70 


0.68 


0.98 


51 


49 


45 


d0 
XJ — > p w 


0.72 


0.39 


0.04 


61 


65 


79 


0.76 


0.37 


0.03 


53 


54 


51 


5° -► iY*°u; 


20.8 


6.97 


1.51 


53 


54 


45 


23.5 


6.58 


1.49 


46 


46 


44 


5+ -»• iY*+u; 


20.8 


6.23 


3.47 


55 


58 


27 


24.3 


6.25 


3.24 


51 


59 


18 


5+ -► 


35.2 


18.8 


1.27 


62 


65 


114 


42.9 


17.9 


1.20 


54 


58 


101 


5° -► iY*V 


33.2 


17.7 


1.20 


62 


65 


113 


40.4 


16.9 


1.13 


54 


58 


101 


B + -► 


0.10 


0.004 


0.67 


59 


1.9 


42 


0.10 


0.004 


0.64 


50 


1.9 


36 


5° -► pV 


0.05 


0.002 


0.32 


59 


1.9 


42 


0.05 


0.002 


0.30 


50 


1.9 


36 


B° ^L0(j) 


0.05 


0.002 


0.32 


59 


1.9 


42 


0.05 


0.002 


0.30 


50 


1.9 


36 
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TABLE XII. B — > branching ratios (in units of 10 6 ) using the BSW form factors in 
models I and II, assuming M H + = 200 GeV, iV e// = 2, 3, oo and tan (3 = 2. 



Channel 


Model I 


SB [%] 


Model II 


SB [%] 




2 


3 


oo 


2 


3 


oo 


2 


3 


oo 


2 


3 


oo 


B° -> p+p" 


17.7 


19.8 


24.4 


0.0 


0.0 


0.0 


17.7 


20.2 


24.3 


-0.1 


-0.1 


-0.1 


5° -> p°p 


0.53 


0.10 


1.21 


0.1 


0.6 


0.1 


0.52 


0.08 


1.20 


-1.8 


-15 


-1.1 




0.90 


0.16 


0.93 


0.6 


1.1 


0.0 


0.86 


0.14 


0.93 


-4.2 


-9.9 


-0.2 


5 + -> p+p U 


13.5 


10.7 


6.04 


0.0 


0.0 


0.0 


13.6 


10.1 


6.04 


2.1 


0.0 


0.0 


B + -> p+w 


14.6 


11.6 


6.50 


-10.9 


-9.6 


-7.4 


14.6 


10.9 


6.50 


-11 


-11 


-7.4 




6.42 


7.13 


8.66 


2.1 


1.9 


1.6 


5.70 


6.25 


7.71 


-9.3 


-10 


-9.6 


s° -> if *°p 


1.70 


1.97 


2.72 


0.3 


0.5 


0.7 


1.47 


1.70 


2.36 


-13 


-16 


-13 




5.44 


6.08 


7.63 


2.8 


2.8 


2.7 


4.90 


5.55 


6.83 


-7.5 


-7.0 


-8.0 


B + -> K*°/9+ 


6.55 


8.13 


11.8 


1.8 


1.9 


2.1 


5.78 


7.51 


10.5 


-10 


-10 


-9.3 




0.37 


0.46 


0.67 


1.8 


1.9 


2.1 


0.33 


0.43 


0.60 


-10 


-10 


-9.1 


^ K K*^ 


0.35 


0.45 


0.65 


1.8 


1.9 


2.1 


0.32 


0.41 


0.57 


-10 


-10 


-9.1 


B U -> p U W 


0.41 


0.23 


0.02 


1.1 


1.0 


0.4 


0.37 


0.21 


0.02 


-10 


-12 


-12 


5° -► K*°u; 


12.3 


4.15 


0.92 


1.3 


1.3 


1.3 


11.1 


4.07 


0.83 


-9.3 


-9.9 


-8.2 


5+ -► K*+u; 


12.6 


3.85 


2.39 


1.6 


1.4 


1.6 


11.3 


3.52 


2.24 


-9.4 


-11 


-4.9 


5+ -► 


20.4 


10.8 


0.65 


1.4 


1.4 


2.2 


18.0 


10.0 


0.54 


-11 


-12 


-16 


B° -► K*V 


19.2 


10.1 


0.61 


1.1 


1.0 


0.5 


17.0 


9.41 


0.51 


-11 


-12 


-16 


5+ -► p+0 


0.1 


0.01 


0.44 


-0.7 


9.7 


2.4 


0.05 


0.005 


0.40 


-11 


11 


-7.5 


5° -► p°0 


0.02 


0.002 


0.21 


-0.7 


9.7 


2.4 


0.02 


0.002 


0.19 


-11 


11 


-7.5 


B° 


0.02 


0.002 


0.21 


-0.7 


9.7 


2.4 


0.02 


0.002 


0.19 


-11 


11 


-7.5 
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FIG. 1. Typical one-loop Feynman diagrams for the quark level decays b — ► (s,d)V* 
(V = 7, Z°, g), with W ± (internal wave lines) and charged-Higgs exchanges (internal dashed lines) 
in the SM and two-Higgs-doublet models. The internal quarks are the upper type quark u, c and t. 



1 



39 



CD 
I 

o 



£ 10 



I 

m 

erf 






100 150 200 250 300 

M (GeV) 






0.0 0.1 0.2 0.3 0.4 0.5 

1 /N eff 



FIG. 2. Branching ratios B(B -> vr+vr^) versus M H + and 1/N eff in the SM and 2HDM's. 
For (a) and (b), we set N e ^ = 3 and M H + = 200GeV, respectively. The four adjacent curves 
are the theoretical predictions in the SM and models I, II and III respectively. The band 
between two dots lines shows the CLEO data with la error: B(B —> 7r + 7r~) = (4.3jr^g) x 1CT 6 . 
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FIG. 3. Branching ratios B{B -► i-C+vr ) versus M H + and l/N e ^ in the SM and 
2HDM's. For (a) and (b), we set N e ^ = 3 and M H + = 200GeV, respectively. The 
dot-dashed, short-dashed, long-dashed and solid curve correspond to the theoretical predic- 
tions in the SM and models I, II and III, respectively. The theoretical uncertainties are 
not shown here. The band between two dots lines shows the CLEO data with 2a errors: 
B{B K+ir°) = (11.6l|;g) x 10~ 6 . 
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FIG. 4. Same as Fig.[|, but for the decay B — > K + ir~. The dots band corresponds to the 
CLEO data with 2a errors: B(B -»■ K+tt") = (17.2±|;|) x 10" 6 . 
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FIG. 5. Same as Fig.|||, but for the decay B — * K°-ir + . The dots band corresponds to the 
CLEO data with 2a errors: B(B -> K°ir + ) = (18.2±|;|) x 10~ 6 . 
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FIG. 6. Same as Fig|| but for the decay B — > K°ir°. The dots band corresponds to the 
CLEO data with la error: B(B K°tt°) = (14.6±|:f) x 10 -6 . 
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FIG. 7. Plots of B{B+ -► if+r/) versus M H+ and l/N e ^ in the SM and 2HDM's. For (a) 
and (b), we set N e ^ = 3 and M H + = 200GeV, respectively. The dot-dashed curve and the 
closely adjacent short-dashed curve refer to the theoretical predictions in the SM and model I; 
while the long-dashed and solid curve correspond to the theoretical predictions in models II and 
III, respectively. The theoretical uncertainties are not shown here. The dots band corresponds 
to the CLEO data with 2a errors: B(B + -► K+r]') = (80l^) x 1(T 6 . 
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FIG. 8. Plots of B{B + -► K+rf) versus M H+ in the SM and model III. The dot-dashed 
line shows the SM prediction with N e ^ = 3. The short-dashed, long-dashed and solid curve 
correspond to model III predictions for N e ^ = 2, 3, oo, respectively. Other theoretical uncer- 
tainties are not shown here. The dots band corresponds to the CLEO data with 2a errors: 
B(B+ K+rf) = (80+|) x l(r 6 . 
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FIG. 9. Same as Fig]^, but for the decay B — > K i]'. The dots band corresponds to the 
CLEO data with 2a errors: B(B° -> K°r]') = (89^) x 1(T 6 . 
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FIG. 10. Same as Fig.[| but for the decay B — > K°ri' . The dots band corresponds to the 
CLEO data with 2a errors: B{B° -> K°rj') = (89+^) x l(r 6 . 
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FIG. 11. B(B -► K* + r/) versus M H + and l/iV e ^ in the SM and 2HDM's. For (a) and 
(b), we set N e ^ = 3 and M H + = 200GeV, respectively. The dot-dashed, long-dashed and 
solid curve show the theoretical predictions in the SM and models II and III, respectively. The 
theoretical uncertainties are not shown here. The dots band corresponds to the CLEO data 
with la error: B(B+ -► K*+rj) = (26.4±^° 8 2 ) x 1(T 6 . 
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FIG. 12. Same as Fig. 11, but for the decay B — > K 77. The dots band corresponds to the 



CLEO data with la error: B(B° -> iT^r/) = (13.8±^) x 10 
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FIG. 13. Same as Fig. 11, but for the decay B — > K* + n . The dots band corresponds 



the CLEO data with la error: B(B 
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FIG. 15. B(B -► versus M H+ and iV 6 ^ in the SM and 2HDM's. For (a) and (b), 

we set N e ff = 3 and M H + = 200GeV, respectively. The upper dots line shows the CLEO upper 
limit: B(B — ► K*°lo) < 19 x 10 -6 . The dot-dashed, long-dashed and solid curve correspond 
to the theoretical prediction in the SM and models II and III, respectively. The theoretical 
uncertainties are not shown here. 
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